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Abstract 


As part of continuing research into the flow about slender bodies of revolution, a con- 
ing motion apparatus for hydrodynamic model testing was built and demonstrated. This is 
the first known use of a rotary balance apparatus for external flow hydrodynamic 
applications. The rotary rig allows for captive model testing with simultaneous roll and yaw 
components, and with non-planar cross-flow effects. Coning characterizes the non-planar 
nature of a general motion that cannot be constructed from contributions due to any planar 
motions. 


Demonstration tests were conducted with a length/diameter = 9.5, body of revolution. 
Force and moment measurements were taken via an intemally mounted, 6-component bal- 
ance. capable of capturing all six force and moment components acting on the model. The 
tests covered coning angles to 20°, spin rates to 200 rpm, and free stream velocities to 30 ft/s. 
Reynolds number based on model length ranged from 4.04x10° to 6.06x10°. Cross-flow 
Reynolds number based on body diameter extended from 1.5x10° to 3.0x10°, covering flow 
regimes near the transition from laminar to turbulent separation. The non-dimensional forces 
and moments generally show a non-linear variation with spin rate for coning angles greater 
than 14°. The measured out-of-plane force was significant, reaching a magnitude of 30% of 
the in-plane force. Experiment data correlated fairly well with numerical predictions (using 
laminar separation criteria). 
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LIST OF NOMENCLATURE 


Body Fixed Axis System 





Z,W,r 
TaN 


eZ. 


G) 


Center of Rototion 


Model Rototion Axis = 





AXES 
Body Fixed: 
26.4 the longitudinal axis, directed from the after to the forward 
end of the model 
ey the transverse axis, directed to starboard 
Od the normal axis, directed from top to bottom 
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Tunnel Fixed: 


° Ap, VR> <p 
*AGs VG. 2G 


~ (06 


c 


me 1, 2 


°K,M,N 


oW 


°B 


me 1,1, 


° CB 
°CG 


the fixed longitudinal axis, colinear with the tunnel longitudi- 
nal centerline. 


the fixed transverse axis, perpendicular to x, in a horizontal 
plane, directed to starboard 


the vertical axis, directed downwards 


COORDINATES AND DISTANCES 


coordinates of the CB relative to body axes 
coordinates of the CG relative to body axes 


coning angle, measured from the x, axis to the x axis 


FORCES AND MOMENTS 


force components relative to body axes, referred to as longi- 
tudinal, side, and normal forces respectively 


moment components relative to body axes, referred to as rol- 
ling, pitching, and yawing moments, respectively, referenced 
to the COR 


weight of the model; W = mg 


model buoyancy force 


INERTIA CHARACTERISTICS 


moments of inertia of the model about the x, y, and z axes, 
respectively 


mass of model 


SPECIFIC POINTS ON BODY 


center of buoyancy of model 


center of mass of the model 
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e COR 


BMC 


Center of rotation, the reference point for all measured and 
calculated forces and moments, located 11.31 inches aft of 
the model bow, on the longitudinal axis 


Balance moment center, .895 inches aft of the COR, on the 
model longitudinal axis 


VELOCITIES 


components along body axes of velocity of origin of body 
axes relative to fluid 


angular velocity components relative to body axes x, y, Z; 
angular velocities of roll, pitch, and yaw 


model rotation vector in tunnel fixed axes; = wu, 
MISCELLANEOUS 


model maximum cross sectional area; area projected onto y-z 
plane 


density of fresh water = 1.9348 slugs/ft’ (t=77°F) 
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A CONING MOTION APPARATUS FOR HYDRODYNAMIC MODEL TESTING 
IN A NON-PLANAR CROSS-FLOW 


Chapter 1 Introduction 


For many years, the primary goal of researchers in the motion dynamics field has been 
to develop the ability to accurately predict the full-scale motions of vehicles. Even today, 
with the powerful computational tools available, reliable motion predictions of vehicles in all 
portions of the maneuvering envelope are not possible. Murray Tobak and Lewis B. Shiff 
pointed out that the difficult problem is to correctly describe the relationship between the 
aerodynamic reactions and the motion variables in the inertial equations of motion of an air- 
craft [1]. This same difficulty applies to the motion of hydrodynamic vehicles as well. For 
slender body shapes with fin appendages, the problem 1s compounded by the complex wake 
structure formed by the maneuvering vehicle. The vorticity shed by the hull and appendages 
creates a wake field that interacts with the velocity distribution over the vehicle’s surface. 
This in tum effects the surface pressure distribution and thus, when integrated over the 
body's surface, the total force on the hull/appendage combination. It 1s this interaction that 
prevents a closed-fomnn analytic solution to the problem. 

To further the understanding of the basic flow field about a slender body, the MIT 
Marine Hydrodynamics Laboratory has been conducting expenmental studies on slender 
body of revolution configurations, both with and without fin appendages. Previous research 
by Coney[2], Kaplan[3], Reed[4], and Shields[5], has focused on experimentally determining 
the nature and strength of the vortical wakes shed from a body of revolution, both with and 


without a single attached fin. Their test program included both force block measurements of 
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fin force and moment as well as laser-doppler velocimeter (LDV) measurements for produc- 
ing cross-flow velocity maps at different stations along the body and for determining vortex 
strengths. 

As both a continuation of the above past research and a step in a new direction, the goal 
of this study is to construct and demonstrate a coning motion apparatus for hydrodynamic 
model testing. The pmnmary purpose of this research 1s to determine the forces and moments 
on the model as a function of the coning rate. The need for the coning motion captive model 


test will be developed from both flow field and mathematical model considerations. 


1.1 Coning Motion 


Coning motion can be described 
as the continuous rolling motion 
of the vehicle longitudinal axis 
about the free-streamn velocity 
vector. To generate such a 
motion in the water tunnel, the 
model 1s fixed to some type of 


support system that can be 





rotated at constant rotation speed 


Vv 
about an axis that 1s parallel to Figure 1. Coning Motion 


the free stream velocity vector of the tunnel. In accomplishing this, the model sees a 
constant attitude with respect to the free stream throughout a rotation cycle. This is a 
Steady coning motion as shown in figure 1. The motion allows for combinations of roll 


(p’) and yaw (r’) velocities simultaneously: p’ =(cos@,)w’, r’ = (sin a)’. 
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In a coning motion study, the experimenter has control over six variables: 


. V: free-stream velocity 
. Ww: body rotation rate 
° @.: coning angle 


. B: sideslip angle 

. : body roll angle 

. d: perpendicular distance of the body reference point from the rotation axis 
For this study, only three of the available variables were used; V, w, and a. The 


sideslip angle, B, and the body roll angle, 6, were set = 0. No sideslip simplified the 
experiment apparatus and simplified the motion for this test. The roll angle is only signifi- 
cant for bodies with attached appendages. The reference point chosen for the experiment 
is the center of rotation, therefore d = 0. 


The choice of the parameter d has a significant umpact on the actual inflow velocity 


to the model. The reference point and the perpendicular distance from the rotation axis for 
this experiment were chosen to minimize inertial force effects, as shown later in section 
2.3.4 and Appendix B. To be completely standard, the reference point should be a function 
of the body geometry (for comparison with other body shapes). The geometric reference 
point is the body center of volume, located on the body longitudinal axis, .8 inches forward 
of the center of rotation. 

The coning motion apparatus has been 1n use by aeronautic researchers since 1926. 
For their purposes, the coning motion apparatus provided test data on aircraft maneuver- 
ability at high angles of incidence and on aircraft steady state spin motion and spin recov- 


ery. Later on, the mathematical model developed by Tobak and Schiff indicated that 
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coning motion is one of the fundamental characteristic motions required for prediction of 
conventional, nonspinning maneuvers. A good historical account of the development of 
the coning motion apparatus and its employment at a number of test facilities around the 


world is given in Reference [6]. 


1.2 Model Testing 


Current hydrodynamic captive model testing comprises planar motion simulations, 
mostly because of the facility limitations. Conventional captive model test types include: 
¢ Straight tow tank or circulating water tunnel at a fixed angle of attack 
¢ Straight tow tank or circulating water tunnel with an oscillator attached to the 
model (1.e., planar motion mechanism) 
¢ Rotating arm 
All of these tests suffer two major deficiencies: 
¢« Cross-flow velocity vectors along the body length lie in a plane 
¢ No roll component of angular velocity (p) 
The significance of these shortcomings can be seen by first looking at the non-planar 


cross-flow effects and then the roll component effects. 


1.2.1 Non-Planar Cross-Flow Effects 


A body of revolution with no attached appendages in a planar cross-flow will 
exhibit four basic flow types, depending on the angle of attack and Reynolds number. 
Shown in figure 2 are the four flows along with the normal and side force experienced 


by an ogive-cylinder at varying angles of attack [7]. 
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The cylinder sheds symmetric vortices with no resultant side force for angles of 
attack up to the onset angle of attack where the flow transitions to steady asymmetric 
vortex flow with substantial side force, and then to wake-like flow with minimal side 


force. 


WORMAL SIDE 
FORCE FORCE 


VORTEX VORTEX FREE FLOW 
FREE FLOW 


SYMMETRIC SYMMETRIC VORTEX 
VORTEX FLOW FLOW SS | | 


a 
= STEADY ASYMMETRIC 

ASY MRE TRIC VORTEX FLOW 

VORTEX FLOW 


WAKE-LIKE Hy 
FLOW WAKE-LIKE FLOW | / | /N 
ANGLE OF ATTACK 


Figure 2. Flow types and force regimes for increasing angle of attack 





on ogive-cylinder. 

Expenments conducted at MIT by Shields[5] on a submerged body of revolution 
at a Reynolds number of approximately 6 x 10° based on model length and at moderate 
angles of attack (up to 14°) demonstrated that the body shed symmetnic vortices for all 
angles of attack. The plot of perturbation velocities for a= 14° shown in figure 3 clearly 
depicts the two symmetnic body vortices. The symmetry here is a direct result of the 


planar cross-flow that the model is exposed to. 
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Some work has been 
accomplished on captive 
models subjected to non- 
planar cross-flows. Visual 
studies done by Tobak, 
Schiff, and Peterson [8] on a 
body of revolution in coning 
motion clearly portrayed the 
asymmetnic vortex field shed 
by the body in the non-planar 


cross-flow. 





Figure 3. Unappended body of revolution, a = 14°. 


More recently, Nielsen Engineering and Research (NEAR) has run a coning 
motion case with their SUBFLO, vortex cloud computer code and predicted similar 
results [9]. The preliminary results are shown in figure 4 for a length/diameter = 9.5 
body of revolution in coning motion with a=20°. The predicted out-of-plane force (i.e., 


side or Y force) was §0% of the normal (Z) force value. 
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Current model test techniques cannot 
capture this effect because of their pla- 
nar flow limitations. Since actual sub- 
mersibles rarely operate in a planar 
cross-flow, the coning motion test 
provides a more realistic test condition 
for obtaining coefficient data. Coning 
charactenzes the non-planar nature of a 
general motion that cannot be con- 
structed from contributions due to any 
planar motions. The presence of the 
non-planar motion is a prerequisite for 
the existence of the coupling terms that 
are sO important in accurately descnbing 
the general motion of a six degree of 
freedom body (e.g., yaw-pitch cou- 


pling). 





Figure 4. L/D = 9.5 hull in coning motion; 


PAVE fol hoy, ee BNO 
1.2.2 Roll Component Effects 


The Taylor Series expansion approach that is described in section 1.3 produces 
many coefficients with a roll component, p’. In the past, these hydrodynamic deriv- 


atives have been ignored, mostly because of the inability to obtain model test data for 
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these derivatives. Three pertinent examples are the linear derivatives, Y’), K’. and 
Nae for which there is little or no experimental results available [10]. eee the linear 
roll damping, is believed to be important on the basis of theory, and has been shown in 
model tests of missiles to be an even function of angle of attack and to vary considerably 
for angles of attack above 5° [10]. 
There are two classes of higher order hydrodynamic derivatives involving roll that 
are of importance: 
¢ Nonlinear variations of roll damping rate with angle of attack: 
eR ar TC oe, 
¢ Nonlinearities associated with cross-flow drag: 
° eee Co. IS ae Mi i Ne 
The nonlinear rol] damping tenns will have a much greater influence on an axisymmet- 
ric body with an attached fin than the body used in this experunent. The cross flow drag 


terms are associated with the nonlinear variation of the cross force with the model angle 


of attack. 


1.3 Hydrodynamic Math Model 


The forces and moments acting on a submerged vehicle are generally non-linear 
functions of the linear and angular displacements, velocities, and accelerations of the 
vehicle and the motion of the control surfaces relative to the fluid [11]. Ideally, the func- 
tional relationship between the forces and moments and the motion and contro! parameters 
would be known and used directly in forming the equations of motion. Unfortunately, the 
functional relationship is generally unknown. Without this relationship, the function is 


expressed as a Taylor Series expansion with respect to the rectilinear and angular velocity 
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components about a chosen condition. For this study, the chosen condition is straight and 
level operation at the instantaneous surge velocity, u(t). Appendix A details a 3™ order 
Taylor Senes expansion and simplification for representative forces and moments. 

For any mathematical representation of a physical process, the validity of the mathe- 
matical model depends on whether the assumed form of the equations of motion ade- 
quately represent the physical hydrodynamics and whether the resulting coefficients are 
realistic. It is not the purpose to this study to validate the use of a Taylor Series expansion 
as a mathematical modeling tool for this case, but rather to present the framework for the 
use of the force and moment data obtained from experiment. The interrelation of the math- 
ematical model, theoretical work, and model test data is best described by figure 5. The 
experimenter s purpose is to provide the experimental results (e.g., non-dimensional forces 
and moments) and the understanding of the dynamics involved to validate the theoretical 


work, thereby improving the comprehension of the actual physical processes at work. 





Theoretical 





Calculations 






Determination 


— | ea 
| Hydrodynamic ‘Characteristic 










~. | | Motions Hydrodynamic 
ode | Information 
___ Sa Model 
Testing 


Figure 5. Mathematical model process. 


The Taylor series expansion is advantageous from the standpoint that it 1s numer- 


ically efficient and provides a convenient structure for correlating data. However, it has 
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the disadvantage in that truncation of the series to a certain order limits the range of 
applicability from the expansion point. Current motion simulators require expansions up 
to 4" order to adequately model yaw-pitch coupling and 5™ order for slender body lift. In 
addition to the algebraic burden of the high-order expansion, for the formulation to be use- 
ful, the coefficients must be accurately determined. 

A simplified expansion for the X force equation, retaining only linear terms is as 
shown: 

XX +X we+X" p’ eX +X" r')+... 
The typical coefficient of a linear term in the expansion takes the form of a partial 


derivative of a force or moment component with respect to a variable evaluated at the orig- 


OX 
x.-(5) 


The form of the hydrodynamic formulation is determined by the coefficients resulting from 


inal condition; for example, 





the expansion. The coefficients determine the characteristic motions that must be evalu- 
ated, either analytically or experimentally. Because of the inability to analytically deter- 
mine the coefficients, model testing, and more recently, full scale ship trials have been 
used for determining not only the magnitude of the coefficients, but also which coefficients 


are of pnmary importance. 
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Chapter 2 Description of Test Apparatus 


All testing reported in this paper was performed in the MIT Marine Hydrodynamic 


Laboratory’s variable pressure water tunnel. 


2.1 Model 


Balance Sleeve 


2.695 = 


| Lf 4.750 ————_ 3.650 Lt 
J — --— Lpmb —————| Dimensions in inches 





Figure 6. Coning Motion Model 


The model used for the experiment has the following characteristics: 


° Length/Diameter 9.5 

. Length 23.50 inches 
. Diameter 2.695 inches 
: Weight 10.45 lbs 

. Buoyancy 3.23 Ibs 

° Jc Parallel Midbody 39.6% 

. Length Forward 4.85 inches 
. Length Aft 10.25 inches 
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The model shape shown in figure 6 is the standard L/D = 9.5, body of revolution 
used for previous studies carried out at MIT (Ref. [3], [4], & [5]), Stevens Institute of 
Technology, Davidson Laboratory, and Nielsen Engineering and Research (Ref. [9]). The 
shell of the model consists of three anodized aluminum pieces which mate to the stainless 
steel sleeve. The sleeve is mounted to the shell of the balance itself. This sleeve is the only 
support for the model shell and allows the balance to sense the forces and moments 
experienced by the model. Void spaces in the assembly were filled with grease to prevent 
water leakage into the model. At the aft end of the model, some clearance from the sting 
was allowed to provide for model deflections. This prevented "shorting out" the balance, 
but also allowed some water to enter the rear cavity. The inertial force contribution from 
this water was calculated as negligible. The next section describes the balance itself and 


its waterproofing. 


2.2 Balance and Sting Assembly 


Center of Rotation Waterproof Sleeve 





Figure 7. Coning Motion Model and Balance Assembly 


Unlike past models, the present body of revolution is equipped with an internally 


mounted, 6-component balance for measuring hydrodynamic loads (figure 7). The balance 
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provided by General Dynamics, Convair Division is a 1000 pound rated balance, calibrated 
to 150 Ibs, and adapted for water testing. The balance was fit inside a waterproof sleeve, 
with the strain gage bridge wires routed out the back through a watertight seal. The bal- 
ance is hollow and allows for wires from the future fairwater force block to pass through 
the balance and out with the balance wires. Great care was taken to ensure that the fine 
wires coming from the balance were adequately waterproofed. Tygon tubing covered the 
wires from the balance to their exit out the water tunnel. 

The balance consists of nine strain gage bridges: 

. 2 Normal force bridges 

“ 2 Side force bridges 

: 2 Roll bridges 

° 3 Axial bridges 
Only six of the bridges are used at one time for producing force and moment readings. 
The preferred set of bridges was determined during the calibration done by Convair. Cali- 


bration of the balance is discussed in Chapter 3. 


2.3 Tunnel Installation 


A considerable portion of this project was spent on the design, selection, building 
and installation of components. During the design phase, emphasis was placed on simplic- 


ity, reliability, and cost. The resulting system 1s shown in figure 8. 
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CONING MOTION APPARATUS 
IN MIT WATER TUNNEL 


Turning Vanes Shaft Encoder 


6.691:1 Belt Drive 
SO Ring Slip Rings 


be-—— Test Section —- 


Y HP 1160 RPM Motor 





Figure 8. Coning Motion Apparatus 


2.3.1 Calculations 


Because of the first time nature of this type of test in the hydrodynamics field, cal- 
culations of the various loads on the system had to be carried out. The calculations were 
integral to the component selection and design process. These estimates included: 

° Model drag and torque at maximum tunnel velocity and rotational speed 

° Model and shaft drive system speed vanation at different rotational speeds 

fora ‘free’ system 

° Shaft loads and system resonances 


. Beanng loads and wear allowances 
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2.3.1.1 Required Torque 


The model assembly drag and torque were calculated by modelling the body as a 
cylinder and the sector as a flat plate. A 2-D strip theory approach was taken, where 
the model and sector were divided into small strips, and the drag calculated for a 


cross-flow velocity equal to wr, for each strip with a drag coefficient based on 


l 2 
D;= cal 5 p(ar,) 


A, is the area of each strip. Torque was then found from summing up the drag contrib- 


Reynolds number: 


utions from each strip: 


Fe DE rD, 


ee | 
Because of the off-center weight effect that the model and sector produced, the 
weight torque had to be added to the hydrodynamic torque. Weight torque was calcu- 
Jated using the same sectioning as in the hydrodynamic case and a density of water for 


the model and a density of .3 Ib/in® (= ,,..,) for the sting and sector: 


Cos y (pV ir, 


as | 


where V, is the stnp volume. The resultant values for the most limiting case, @, = 


20°, v = 30 ft/s, w = 200 rpm, were t, = 14.9 ft-lbs, t, = 5.9 ft Ibs. Finally, the mini- 
mum required HP was found to be: 


(T, + 4,,)@ 


= .79 HP 
550 ? 


HP 
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Later testing proved two problems with this result. The first was the calculation 
was Carried out for the preliminary sector configuration, which was more streamlined 
than the final design. Second, the effect of the high free stream velocity (when com- 
pared to the wr contribution) was neglected. The addition of the V term caused the 
sector to act as an inefficient foil at a varying angle of attack along its span (tunnel 
radius). Later attempts to model this correctly produced a very flat torque characteris- 
tic curve and also proved wrong. The effect of this error was the limitation of approx- 
imately 135 rpm for the model at 20°, 200 rpm, and 30 ft/s. 200 rpm could only be 
achieved for the model at 20° when the free stream velocity was reduced to 
approximately 10 ft/s. This problem was later fixed by placing a wing attachment on 
the end of the sector. The attachment was designed for a 2 to 3 degree a at 200rpm 
and 30ft/s. The wing was a large success, providing enough lift to autorotate the 
model at 175 rpm (for V = 30 ft/s)! With the motor assisting, the model easily 


achieved 200 rpm. 


2.3.1.2 Speed Regulation 


For rotary balance experiments, it 1S umportant to maintain a constant rotation 
rate over acycle of data taking. To check the worst case condition, a simple calcula- 
tion was done for a freely spinning system (i.e., no speed regulation). The model and 
sector were lumped together as an off-center weight of 26 pounds at a radius of 
gyration equal to 4.24 inches. The main sheave polar moment of inertia was calcu- 
lated and later obtained from the manufacturer. All other polar inertia contributions 
were neglected on the assumption that they were small. 


For this lumped system the total J was: 
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ln ae (J model a5 Tee) = 4366in’ an lbs 


The following relations were used to calculate the rpm variation: 





_ d(H.) = 
ae ay where H,= ) myr;@_ ,the angular momentum for this case. 
sm) 
nee See ee — 
Substituting: eS De, mr,d = Jw 
ia) 


Solving for the angular acceleration w: 
ee 
Clee 7 =A sin(@,') where, , = 20.9rad/s 
Finally, the angular speed variation was found from integrating this expression: 


A 
w= | xt = aes =P (Dp 


oO 


Substituting in the numbers for 200 rpm (20.9 rad/s) and 1,. = 10 ft-lbs: 


w= -—.507 cos(20.9r) + 20.9 


This gave a worse case vanation of = 4.8 rpm over | revolution. The advertised speed 
regulation for the motor would give a max variance of = 2.2 rpm over 1 cycle (1%). 
To avoid all of this, the main sheave was counter weighted during the inertia and data 


runs. 


2.3.1.3 System Torsional Resonance 


Unknown loads on the model-sector combination could result in a resonant 


vibration in the system. To investigate this, the torsional resonance frequency was 
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calculated for varying J’s and shaft "stiffness" constants (K,). The system was mod- 


eled as atwo disk, torsional vibratory system, shown in figure 9. The model/sector 


combination accounted for J,, the main sheave, J,. 





Figure 9. Torsional Vibration Model 
The equation for the resonant frequency was derived starting with Lagrange’s 


equations with no forcing function: 


L Le : 
Equation of motion: “(%|- (= |= O ,where L=T -V 


dt 


f= (: \7,6; -- (= 7.63 (Kinetic co-energy) 
2 


] : 
V= ke —6,) (Potential energy) 
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Assuming an undamped system, the angles take the harmonic form: 
Q,=A,e"™ 
6,=A,e"" 


The resulting equation to solve is: 


eee Sls 
ger ER 


The characteristic equation is: det[B] =0. This leads to the desired relation for : 


Rearranging: 





; ee 


Ro Gs 4 
—— K == = 
1) Tid, | , Where K,=g2 39) (d,-d, ) 


d, = outer shaft diameter, d= inner shaft diameter, /= shaft length, G = 12 x 10° psi 


for steel, and g = 386 in/s’. For the final system, the torsional resonance frequency = 
50 Hz (3000 rpm), which is much higher than the filter cutoff frequency of 3 Hz (180 


rpm) and the maximum shaft rotation rate (200 rpm). 


2.3.2 Modifications 


The MIT variable pressure water tunnel required significant modifications before 
the study could be undertaken. The primary modifications necessary were: 

. Boring and sleeving the turing vanes to allow for the shaft installation 

2 Modifying the forward nacelle to accept "DU" type bushings for supporting 


the shaft, model, and sector in the test section 
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Aligning and welding the bearing housing mount plate on the tunnel exterior 


Installing a motor and drive system for rotating the model and shaft 


Special consideration was given to the rigidity of the system since unknown reso- 


nances in the model-sting-drive system would have a disastrous effect on the data. The 


forward nacelle is nmgidly mounted to the tunnel as is the exterior bearing housing 


assembly. The precautions proved to be worthwhile; the apparatus showed only small 


amplitude vibrations up to maximum system capabilities. 


2.3.3 System Characteristics 


The basic system characteristics are: 


Model angle of attack 0° to 20° in 2° increments 

Rotations speed 0 to 200 RPM, both directions 

Tunnel water velocity 0 to 30 ft/sec 

Rotary speed regulation to + 1% for 95% change in load 

Drive motor: GE 2 HP DC shunt wound, 1200 rpm, variable speed 
Drive system: 2-belt drive with drive ratios from 5.59:1 to 8.2:1 


50 channel (rings) capability 


Model rotation rate was derived for similitude in non-dimensional roll, p’, with 


full scale vessels. Using the conventions in Appendix A, p’ = 1.37 for 200 rpm at U = 


30 ft/sec. The coning angle limit of 20 degrees is based on tunnel test section dimen- 


sions. 
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2.3.4 Test Section Installation 


There are three contributions to the forces and moments measured by the balance 
in the model (see figure 10). The first is from the inertial forces and moments of the 
model itself, which vary with the model attitude and rotation speed. The second 1s from 
gravity (model weight) and buoyancy, which will vary cyclically over a revolution. 
Finally, there is the contribution from the hydrodynamic loads. The desire to reduce or 
eliminate all but the desired hydrodynamic forces had a great effect on the installation 


design. 


Gravit y: Periodic 
Inertia: Constont 
Hydrodynamic: Constont 


iz, Weight (or Buoyancy) 


Force 


or 
Moment Hydrodynamic 


Inertial 





Figure 10. Forces measured by the balance for constant rotation rate 


A close-up of the model in the test section 1s shown in figure 11. The model cen- 


ter of rotation was chosen such that the model center of gravity would be as close to 
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Center of Rotation 


Mode) Rotation Axe -_ 3 5 
Sector Support oe 


Test Section Wall 


Dimensions in Inches 





Figure 1]. Coning Motion Apparatus in the Test Section 


the model rotation axis as possible. This reduces the inertial force contribution. Addi- 
tionally, reducing the differences in the model moments of inertia also decreases the 
body inertial moments (see Appendix B). As it turned out, the inertial forces were very 


small compared to the hydrodynamic loads. 


2.4 Data Acquisition and Reduction System 


The data were taken with the stand-alone, microcomputer-based, data acquisition and 
reduction system shown in figure 12. The components wull be descnbed by their func- 


tional grouping. 
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Figure 12. Data Acquisition and Reduction System 


2.4.1 Instrument Group 


The instruments consisted of the 6-component balance, tunnel differential pressure 
(DP) cell, and shaft encoder. Their basic functions were as follows: 
. Balance: 9-element (Strain-gauge bndge), internal balance used for measur- 
ing the forces and moments expenenced by the model 
- DP Cell: Used for measuring test section velocity 
. Shaft Encoder: 12 bit, natural binary encoder (BEI model C-14) used for 
model rpm, data tngger, and shaft position for static tests. The encoder out- 


puts a 0 to +5 volts square wave signal at frequencies from 1/revolution to 
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2048/revolution (2''). The 1/rev and 32/rev signals were used for shaft 
position, the 32/rev for tniggering, the 512/rev for model rpm (computer 


input), and the 1024/rev for model rpm (counter input). 


2.4.2 Amplifier/Filter Group 


Signal amplification and filtering was accomplished by two sets of instruments: 

° 3-B Series (Type 3B-16): A nine module block of signal conditioners and 
filters that filter and amplify the balance outputs to a standard + 5 Volt ana- 
log signal. The modules had an upper cutoff frequency of 3Hz for filtering 
out high-frequency noise (a characteristic that caused considerable trouble 
later). 

> Daytronix 9000 Series: Two modules were used: the strain gauge condi- 
tioner (9170) for boosting the DP output to a+ 5 Volts range, and the fre- 
quency to voltage module (9140) for converting the encoder 512/cycle 


output to aQ to +5 Volt analog signal for model rpm. 


2.4.3 Computer Group 


The heart of the system was the IBM PC/XT personal computer. The PC had a 
MetraByte Dash 16F analog to digital (A/D) board installed for acquiring and converting 
the balance voltages to counts. The Dash 16F also read in outputs from the tunnel dif- 
ferential pressure (DP) cel] and the shaft encoder. The analog portion of the board was 
set for a resolution of 2.44 millivolts (409.6 counts/volt). The board had a conversion 
tume of approximately .0417 mullisec./channel using direct memory access (DMA). This 


was a very important point which will be discussed in chapter 3. 
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The Dash 16F also was capable of reading in a digital signal through its digital 
input ports. This feature enabled the encoder to act as a trigger, using software to deci- 
pher the digital word read in from the port. 

The PC/XT controlled the data taking process through software. Standard Project 
Athena laboratory routines were used for controlling the analog and digital sampling. 
The routines were incorporated into user-developed FORTRAN code. The major draw- 
back to using the FORTRAN routines was the increased time required to obtain an ana- 
log sample. The experimentally determined sample time was approximately .052 
msec/channel, an increase of .0103 msec/channel over the calculated rate. 

Finally, the PC stored and processed the averaged component (DC) of the signal 
and on-line presented the force and moment measurements. The time-varying signal 
was stored on magnetic media for later off-line processing. 

Not connected to the computer, but used for setting the rotation rate was the HP 
(Hewlett-Packard) counter shown in the figure. The counter read directly off the shaft 


encoder and gave an accurate counts/sec output, which was then converted to RPM. 


2.4.4 Slip Rings 


The balance signals (millivolts) were taken from the rotating model reference 
frame to the fixed data system frame via a 50 ring slip ring assembly. The slip ring 
assembly (Airflyte model DAY-491-50) used was specifically designed for strain gauge 
measurements. As the experiment showed, the rings were very “quiet”, passing little 


electrical noise to the amplifier/filters. 
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Chapter 3 Test Procedure and Tare Measurements 


3.1 Balance Calibration and Check-out 


The 6-component balance calibration was accomplished at Convair, prior to its 
installation in the model. A precise calibrating body was fit over the balance and then the 
assembly put in a calibration fixture. The fixture allowed for the accurate placement of 
loads and moments on the body, with corrections for deflection being taken into account. 
The calibration results were placed in a two inch thick binder containing information nec- 


essary for using the balance. The information included: 


. R-cal readings for all 9 channels (shunt resistor equivalent) 

. Balance constants for the different combinations of roll moment and axial force 
gauges 

° Plots of the response of the balance vs applied load (for applied axial, normal, 


and side force and roll moment) 

The R-cal readings are merely voltage (counts) outputs for the balance channel being 
tested, with and without the shunt resistor switched into the circuit. The A voltage (counts) 
allows a second site to establish calibration ratios between the counts read on their A/D 
system and the equivalent readings on the Convair test bench. In this experiment, the ratio 


Was Set aS: 


: R-Cal caiveir 
Ratio= yey can = 2.45 


R-calygr 
The ratio vaned by channel, however, the vanance was small. The ratio calculated above 
shows that the Convair A/D system had approximately a factor of 2.45 better resolution 


than the system used in this study. 


Paces> 





The balance constants were not used because of the simplified linear algorithm cho- 
sen for converting the balance readings to forces and moments. Several methods were 
available for doing this conversion, but the linear method was chosen for its simplicity and 


speed. Other methods available are described in References [12] and [13]. 


3.1.1 Data Reduction Model 


Stevens Institute of Technology, Davidson Laboratories, developed the linear, 6" 
order method unplemented in this study. In developing this linear model, Stevens had to 
fit the Convair calibration data and then compare the linear model performance against 
the more complex Convar algorithm. 

The first step is to compute the 6x6 calibration matrix. The goal is to obtain a 
matrix of coefficients that accounts for the umperfections 1n the balance construction, 
which appears in the fomn of cross-talk, and yields the correct force vector. In other 
words, for a sunple single. in-plane force, such as a normal force applied at the reference 
center, the detected force by the balance will consist of contributions from all 6 chan- 


nels: 


Ze= aie 
im] 


where a, are the calibration coefficients and r; are the balance readings from the six 


channels. For a generalized load application, the matrix equation is: 


L=[A]R 
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where L is the load vector, [A] 1s the desired coefficient matrix, and R is the correspond- 


ing balance reading vector. During the calibration done at Convair, several load vectors 
were applied while recording the associated balance reading vectors. This resulted in 


the following, overdetermined system: 


Z, Z, Qi, Gyr yy Ay ys Aye{| NI, N1, 
M, M,, G3, Ax, Gz; Ax, Ags Ayg|| N2, N2, 
Vy Dy, | [31 32 433 G3q 35 36 DG Cnr 
N, Np , Bg, Ayn Ag, Agy Ags Age} | Y2, ris 
| K, “i lid Asx G53 Ase Ass As6 ji R2,,| 
Xx a Gp, (Ge Ag, Gun (Aces |), Al! 


where the subscripts refer to the loading condition. This equation can be rewritten as: 
[R]' [A}' =[L]’ 
The solution to this system is found by utilizing a least squares approach. This is equiv- 
alent to minimizing the Euclidean norm of 
[R]" [A}’ -{L]' 
Davidson Lab has a standard program which, given a fule containing the loads and 
results, carries out both the linear least squares fit and the matrix inversion to get the 
final coefficient matrix, [A] [12]. This matrix was then used to convert the balance 


readings obtained from the experiment into the force and moment data presented. 


3.1.2 Balance Operation 


The next logical step was to verify the correct operation of the balance. Successful 
completion of this portion of the experiment would give: 


° A check of the data reduction method 
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‘ The balance axis convention 
. A check of the data taking system 


° Model weight and center of gravity 


3.1.2.1 Weight Calibration and Axis Convention 


The first experiment min was the application of known weights on the balance. 
Stevens Institute provided a collar attachment that fit around the stainless steel sleeve 
and allowed weight placement with the model set at a. = 0 degrees and roll = 0 
degrees. Weights were systematically placed on the collar and readings taken using 
the static FORTRAN code developed for these runs. Applied loads varied from .4 to 
12 pounds. The resultant Z force was approximately 1.8 times the applied load. After 
much searching, the error was found in the coefficient matrix associated with refer- 
ence to the center of rotation. The transpose of the coefficient matrix had been used, 
causing the Z force error. The error was corrected, allowing the next step to proceed. 

The sign convention for the axis system was checked by furst setting the model 
to 0°, 0° and taking a measurement. This measurement became the "zero" file and 
would be subtracted from the subsequent reading. Next, a known (in sign only, not 
magnitude) force or moment was applied to the model and a set of readings taken. 
The difference between the second reading and the zero file produced the read sign of 


the force (or moment). The test showed the following axis convention: 


° X: positive out model bow 

° Y: positive out model port side 

° Z: positive down 

. K: positive for model starboard side up 
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. M: positive for model bow down 

° N: positive for model bow to port 
(Bold typeface indicates non-standard sign convention.) To avoid further confusion, 
the axis convention above was adopted for the test readings only. The test measure- 
ments would be presented as follows: 

° Dimensional forces and moments: above, non-standard convention 

° Non-Dimensional forces and moments: SNAME axis convention as pres- 


ented in the list of nomenclature. 


3.1.2.2 Model Weight and CG 


With the axis convention determined, the model weight and center of gravity 
were established. Two types of experiments were run, first, an inclining experiment, 
second, a simple set of roll experiments. In the first experiment, the model is set at 0° 
(pitch or &.), O° (roll), a set of zeroes taken, and then inclined 2 degrees at a time (in 
pitch) with readings taken at each point. The weight is found from the relation 


X =W sin(a.) 
The X force is plotted against sin(a,.). The slope of the plot is the model weight. The 


slope gave the model weight as = 10.3 lbs. Though the data plotted very linearly, the 
maximum X force measured was = .35 pounds. More emphasis was placed on the 
results of the following roll tests. 

The next experiment consisted of a set of roll runs. First, the model was placed 
at either 0°, 0° or 0°, 90 (degrees roll angle). A set of "zeroes" were taken. Then the 
model was rolled 180 degrees and again sampled. The resulting readings gave the 


following results: 
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AK 
ge” 2) 


+ 000, 180 Test: Z = 2(Weight), X,, = =m, Y, 


AN AK 
° 090, 270 Test: Y = -2(W), X., = = La= aa 
Giving numbers: 
° W = 10.45 pounds 
. X,, = -4 inches forward of center of rotation 


> Nae a 

As a measure of some confidence, the test data was re-reduced using the coefficient 
matrix referenced to the balance moment center (BMC). The results were nearly iden- 
tical. Additionally. the different components of the balance gave consistent results 


(1.e., Y force detennined W = Z force determined W). 


3.1.2.3 Temperature Effects 


During this check-out phase, the temperature sensitivity of the balance was dis- 
covered. Because the balance is a large mass of metal, the temperature of the balance 
took some fixed time to stabilize. Even the smal] heat generated by energizing the 
gauges affected the balance outputs (significantly). A quick experiment was run 
where, with the model in air, the balance was energized and readings taken every few 
minutes. The first 15 minutes produced drifts of 2.5 Ibs Y, 1.1 Ibs Z, 1.9 in-Ibs K, and 
1 in-lb M. Over the next 30 minutes, the drift fell off considerably. 

To counter this "warm-up’ effect, the strain gauges were left energized during 
the entire experiment. The temperature sensitivity caused problems in the full-up data 


runs later because of the tunnel water temperature mse over mun time. 
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3.2 Test Procedure 


3.2.1 Test Logic 


A coning motion test produces forces and moments from multiple sources. Three 


components discussed in Chapter 2 are weight (or weight - buoyancy), inertia and 


hydrodynamic. An additional "force" is the bridge offsets or zeroes. This fourth com- 


ponent is a function of temperature. The only force (or moment) of interest is the hydro- 


dynamic force. To this end, a test procedure had to be developed that accounted for 


each component and left the desired result. For this study, the desired result was both 


the time varying (or AC component) and the steady state (or DC component) hydrody- 


namic force. The basic test procedure was (starting with the model at the desired @. ): 


A. 


1) 


2} 


Model in air 


Take a set of readings at set angular positions around a 360° rotation. Mea- 
surements done every 11.25° starting at O° and rotating in the positive (stan- 
dard convention) direction (by hand) to 348.75°(32 sampling positions). 
Store the raw counts for all 9 channels (100 sample sweeps/angular 
position, summed to | point), both for each point (AC data) and the summa- 
tion of all 32 points (1 revolution)/channel (DC data). Result: weight and 
offsets raw counts for that @, . 

Conduct "wind-off' runs at same rotation rate (@) and direction as "full-up" 
runs. Sample all 9 channels (1 sample sweep/angular position) and store 
each raw counts data point. Sum raw counts for each channel over # of 
revolutions (usually 10) and store for processing with DC water nun results. 


Result: weight, offsets, and inertia raw counts for that @, and . 
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3) Subtract DC raw counts measurements in test A.1 from the DC raw counts 
measurements in test A.2. Process the results for only the preferred 6 chan- 


nels. Result: inertia force for that o, and @. 
4) Repeat A.2 through A.4 for each w, A.1 through A.4 for each @, . 


Note: Though not processed immediately, the weight and inertia vs angular posi- 
tion raw counts are stored and can be processed later for the inertial effect 
over a rotation cycle. 

Model in water 

1) After a “soak” period, usually overnight, repeat test procedure A.1, with the 
model in water and tunnel free stream velocity = 0. Result: weight - buoy- 
ancy, and offsets. 

2) Establish tunnel test conditions (1.e., tunnel free stream velocity and model 
rotation rate). Conduct "full-up" run, recording |1 channels of data (9 bal- 
ance, | DP, 1 rpm) at each angular position and storing the raw counts. 
Sum raw counts for each channel over # revolutions for averaged (DC) 
component. Result: weight - buoyancy, offsets, inertia, and hydrodynamic 
raw counts. 

3) Subtract DC raw counts measurements in test B.] from DC raw counts 
measurements in test B.2 (9 balance channels only). Process averaged data 
for only preferred 6 channels, DP cell, and rpm data. Result: inertia and 


hydrodynamic averaged force (over # rotation cycles). 


4) Subtract inertia force measurements corresponding to that ® and a. from 


the B.3 result. Result: Hydrodynamic force for those test conditions. 


5) Repeat B.2 through B.4 for each w, B.1 through B.4 for each a, . 
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Note: Again, the tume varying data are all stored as raw counts files for future 
processing. 
The motivation behind the many decisions made in preparing this procedure is 


discussed in the following sections. 


3.2.2 Code Development 


At the core of the data taking was the FORTRAN software developed by the 
author (in conjunction with the routines developed by Glenn McKee of Stevens Institute, 
Davidson Labs). The form of the test procedure and its limitations dictated the form of 


the code. The critical issues considered were: 


° Timing 
° System capabilities 
° Desired outputs (i.e., time varying and averaged forces) 


3.2.2.1 Timing 


The basic timing parameters considered were: 


. A/D sample and conversion time for }] channels 

° Model rotation rate 

° Code loop execution time 

. # of divisions over a revolution 

° Maximum allowed change in gravity vector over a sample sweep 


The MetraByte Dash 16F specifications give an estimated A/D conversion time 


with DMA (Direct Memory Access) transfer of .04167 msec/conversion. For 11 
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channels, this works out to = .46 msec. The experimentally determined A/D conver- 
sion and transfer tume = .052 msec/conversion giving a total = .57 msec for 11 chan- 
nels. The importance of this is seen by looking at the maximum allowed change in 
gravity over a sample sweep. 

Using a small angle approximation for the change in rotation angle over a sam- 
ple sweep, the following relation was derived: 


AG = A®@ = .Olrad (.57 deg) 


for amaximum allowed change of 1%. At the maximum rotation speed of 200 rpm, 
this works out to .4775 msec. The data taking system used simply could not satisfy 
the 1% criteria (at 200 rpm). The actual AG for this study was = 1.2% (.68°). 

The time to execute one data taking loop lead directly to the # of increments that 
could be sampled over | revolution. Use of the shaft encoder as both a trigger and a 
position indicator (locating the 0° position) required several logic statements in FOR- 
TRAN. The logic allowed the code to read the digital word on the D/D input ports 
and determine, |) when to set 0°, and 2) when to tngger a sample sweep. Embedded 
in the logic loop was the Project Athena Laboratory sampling routine for reading in 
and storing the A/D data. By running several loop timing versions of the data taking 
code, the time to execute the loop was found to be = 5.2 msec for an 11 channel 
sweep. Again looking at 200 rpm, the maximum number of increments in one revolu- 
tion that would allow the loop to execute and retum for the next trigger pulse was 32. 
This allowed 9.375 msec for loop execution and retum. 32 data points per rotation 
were felt to be adequate resolution for the purposes of the experiment. 


The data taking code used in this study is included in Appendix D. 
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3.2.2.2 Desired Results 


During the initial code development, the thought was to work directly with the 
individual raw counts data points for subtracting tares and converting to forces. The 
averaged force values would then be computed from these results. Unfortunately, 
upon actual data taking, a mysterious phase shift and amplitude drop off were noticed 
for all rotational speeds. The offenders in this case were the 3B-16 Amplifier/Futers 
supplied with the balance. This should have come as no surprise, especially when the 
3 Hz upper cutoff frequency was known a priori. The phase shift can be seen graphi- 
cally in figure 13. Fortunately, only the time-varying or AC component was affected, 


passing the zero Hz or DC component unchanged. 
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Figure 13. Phase Shift in Raw Count Data 
This discovery initiated a flurry of code rewntte activity and caused a change in 
philosophy to working only with the averaged counts. To maintain the ability to 


work with the time varying data off-line, and to quantify (at least for one channel) the 
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amplitude reduction and phase shift, a bode plot was made. The phase shift is shown 


in figure 14, the amplitude plot in figure 15. The best fit line for each was recorded 


and stored for later use. 
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Figure 14. Phase Shift vs Frequency 


Page 46 





Ampirtude Rotio vs Frequency, Chonne! 0 
Sth Order Fit 


O 
pe 
Oo 
ar 
w 
o 
= 
os 
or 
5 
< 


0.1 — 
Frequency (Hz) 





Figure 15. Amplitude Ratio vs Frequency 


3.3 Tare Measurements 


The test procedure delineated in section 3.2.1 accounts for the (W-B) and inertia 
effects. Though it is possible to account for these forces analytically, they were eliminated 
in a more straightforward manner by subtracting averaged (DC) measurements. The tare 
tests were run in air since, compared to water, the air acts essentially as a vacuum and 


doesn’t influence the results. The general idea for the conduct of the tare measurements 


came from Reference [14]. 
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3.3.1 Static 


The static (W = Q) tares were very simple. The model (in water) was set at the 


desired a, the roll angle set to 0, and a set of readings taken. This "zero" file accounts 

for the (W-B) and offsets for that a. The tunnel test conditions were set and another set 
of readings taken. Subtraction of the two files left the desired static hydrodynamic load. 
A more thorough approach would have been to repeat this at several rotational positions 


and average the results. Because of time constraints, this was not done. 


3.3.2 Rotational 


The inertial loads can be calculated analytically, given the model mass, moments 
of inertia about the principle axes, and the location of the model CG. Additionally, these 
values must be assumed to not vary with rotation rate. If this course is pursued, then the 
following relations, derived in Appendix B, result: 


¢ For the inemia forces (in SNAME axis convention): 


re? 
X, = mW xX sin’ a, 


Z =-mM@ X¢ Sin &, cos &, 
M =-w'(/, —1)(sin a, cos a,) 
Ka N10 


¢ For the weight and buoyancy forces: 
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Xw_, =—(W -—B)sina, cos wr 
Yy_, =(W -B)sinor 
Zy 3 =(W —B)cosa, cos wr 
Ky. =W[y, cos @, cos @f — Z, sin wr] 
M,,.=—-W cos @r[z, sina, +%,_ cos &] 
Ny =W [x, sin @r + vy, sina, cos Wr] 
The same moment equations result for buoyancy by substituting -B for W and Xp, Vp, Zp 


fOr XG. Vg: Z¢ - 

There are several problems with using these relations directly. Furst, the 
model/sting combination deflects during rotation, causing the actual inertial loads to 
vary. Second, use of the weight and buoyancy relations would involve determining the 
offsets separately, rather than the more direct method of lumping them with the weight 
and buoyancy effects. Finally. the relations depend on accurately knowing the model 
mass moments Of inertia, the CG, and mass. Errors in these quantities would reflect 
directly in the calculated forces and moments. 

For this experiment, the deflection should be small due to the relatively small iner- 
tial forces at the maximum @, and @. As it tumed out, the design of the model/sting 
minimized the inertia force contribution. The calculated and actual measured inertia 
forces for a = -20° and -14° are shown in the following figures. Figures 16 and 17 are 
for the -20° setting, figures 18 and 19 for the -14° angle. The moment measurements 


have been corrected for sign to conform the standard convention. 
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Figure 18. X and Z Inertia Forces, 14° 
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Because calculation of I, and I, would be approximate at best, the actual data was 
fit with a second order curve, and the result used to empirically determine the two terms. 
This method worked quite well, with only a 2.9% difference between the 14° and 20° 


curve fit determined values. The numbers used in calculating the forces and moments 


wele: 
‘ Coie -.3706 slugs-in’ 
° m: 3248 slugs 
‘ XG: .40 inches 


The 14° X, Z, M, and the 20° M results seem to match very well with the calcu- 
lated trends. The 20° X and Z force measurements depart significantly from the calcu- 
Jated trends. Much scatter can be seen in the experiment results. Three factors could 


explain the results: 


- The measured forces are very small, much less than 1% of full scale calibra- 
tion load. 

° The temperature effects were not well compensated for. 

° X¢ and m are off from their actual values. 


The third bullet helps explain the excellent correlation for the M plots and the poorer 
correlation for the force plots. Also, the moment correlation might be better because the 
moment values arise from subtraction of the two normal force balance outputs, possibly 
cancelling any errors. The Z force readings are the result of the addition of the two 
norma! force balance readings, possibly compounding the error. 

Despite the force discrepancies, the actual inertia force is small compared to the 
hydrodynamic force, ranging from 1.7% to 4% of the measured Z hydrodynamic force. 


The only significant inertia effect, M, is well predicted. 
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One point that has been overlooked is the deflection of the balance due to the 
hydrodynamic force. No attempt was made to correct for this effect, and is left as a 


refinement for future experiments using this apparatus. 


3.4 Test Matrix 


Before any rotational testing could be done, the test matrix had to be determined. 


The matrix evolved from satisfying the goals: 


° Widest possible coverage of non-dimensional spin rate, w’. 
° Maximum hydrodynamic force. 

. Minimum Reynolds effects. 

° Maximum # of @ increments In a test mun (i.e., O - 200 rpm). 
’ Maximum @, coverage (1.e., 0 to 20°). 

, Maximum p’ andr’ combination. 


Because of time constraints, only three @, settings were tested; -8°, -14°, and -20°. For 


each a, the tests were run in increments of 12.5 rpm, from 0 to 200 rpm (where possible), 
in one rotation direction, and at 25 ft/s free stream velocity. The variances on this proce- 
dure were: for -20°, do tests for 20, 25, and 30 ft/s free stream velocity to investigate 
Reynolds effects, for -14°, do runs in both rotation directions to get function even/odd 
characteristics. 

The 25 ft/s base velocity was chosen as a compromise between hydrodynamic force 
and range of non-dimensional spin rate. This resulted in a minimum Reynolds number in 
the range of 5x10°, well above the agreed upon lower cutoff of 3x10° (for Reynolds num- 


ber dependencies). The resulting p’ andr’ ranges were: p’: 0 - 1.93, r’: 0 - .70. 
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Chapter 4 Test Results 


"Full-up” test mins were conducted using the rotary rig and standard L/D=9.5 model 
described previously. The test procedure delineated in section 3.2.1 was followed. Normal, 
side, and axial forces, rolling, pitching, and yawing moments were measured with the model 
in coning motion at coning rates up to 200 RPM. The majority of the runs were performed at 
25ft/s (Reynolds Number based on model length of 5x10°). Coning angle, @,, varied from 
-8° t0 -20°. The measured forces are non-dimensionalized and presented as functions of the 
non-dimensional spin rate, @’. The axial force and roll moment are not presented as they 
were determined not to be particularly important for this experiment. All non-dimensional 
parameters are presented in the standard SNAME, body-fixed, axis convention [15]. 

Figures 20 and 21 show the non-dimensional force and moment variation with both 
spin rate and Reynolds number for a, = -20°. Figures 22 and 23 show the non-dimensional 
force and moment variation with both positive and negative spin rates for a, = -14°. Figures 
24 and 25 present the force and moment vanation for a, = -8°. Finally, figures 26 and 27 
show the variation of the non-dimensional side force and yawing moment, respectively, with 


coning angle. A tabular listing of the forces and moments 1s in Appendix C. 
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Figure 20. Y’ and Z’, a, = -20° 
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Figure 21. M’andN’,a, = -20° 
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Figure 22. Y’andZ’,a, =-14° 
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Figure 23. M’andN’,a, =-14° 
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Figure 24. Y’ andZ’,a, = -8° 
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Figure 25. M’andN’,a, = -8° 
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Figure 26. Y’ vs-a, 


The 20° force data demonstrate a nearly linear dependence on coning rate for the nor- 
mal force and a very non-linear dependence for the side force. Of note is the significance of 


the out-of-plane force (Y’), which reaches 30% of the in-plane force. Both the pitching and 
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Figure 27. N’ vs-, 


yawing moments show a non-linear dependence on w’. Reynolds number effects seem negli- 
gible for the forces, and to have a slight influence on the moments. For the yawing moment, 
the lower Reynolds number data shows a tendency to go non-linear sooner than the 25 and 30 


ft/s data. The difference in the pitching moment 1s not readily explained. 
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The normal and side forces have the even-odd characteristics anticipated. Because of 
symmetry, the lateral force and moment, Y’ and N’, should be odd functions of the coning 
rate parameter, while the longitudinal force and moment, Z’ and M’, should be even func- 
tions of the coning rate. The 14° data clearly shows the even/odd relationship for the forces 
and moments. The slight asymmetry in the data with respect to positive and negative spin 
rate could result from misalignment of the model/sting combination with the free stream 
velocity. 

The 8° force data exhibited simular trends as the 14° and 20° results. The pitching and 
yawing moments, however, appear almost constant over the range of spin rates investigated. 
The slight deviations for the zero spin rate results are most likely test anomalies, and not true 
pictures of the trends. The N’ static result should be zero, as a linear interpolation through 
the remaining data points suggest. 

The moments and the normal force display simular trends for the 14° and 20° results. 
The side force shows a different character between the &, settings, especially when the &° 
data is included. As figure 26 shows, the side force 1s a very non-linear function of coning 
angle, changing signs at approximately 14°. The relationship could be quadratic, but with 
only three data points per line, the exact form is impossible to determine. The yawing 
moment also shows a non-linear characteristic with coning angle, becoming nearly linear for 
the high spin rates. Obviously, more data points are necessary before any strong functional 


dependencies can be determined. 
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Chapter 5 Data Analysis 


The purpose of this chapter is to validate the steady-state data presented in Chapter 4. 
Two main paths are pursued: 
. Comparison with numerical results. 


° Frequency analysis of raw balance data and encoder determined RPM data. 


5.1 Comparison with Numerical Results 


A vortex cloud numenical algorithm, developed by Nielsen Engineering and 
Research (NEAR), was used for producing the calculated results for comparison. A 
description of the computer code SUBFLO is found in reference [16]. The method basi- 
cally models the body and major physical features of the flow field with singularity distnb- 
utions. Mutual interactions of the body shed vorticity are considered in the prediction of 
the induced forces and moments [9]. The test cases shown in the following table were 
provided to and run by NEAR for producing the predicted forces and moments. The 
NEAR predictions were made for both no separation and laminar separation. Laminar sep- 
aration was used to avoid the uncertainty of transition effects. 


In all cases, except for the 20 ft/s, 125 rpm, N’, result, the experimentally determined 


force or moment fell in the range bounded by the laminar and no separation cases. Addi- 
tionally, the laminar separation cases more accurately predicted the experimental result for 
all runs. For this study, the cross-flow Reynolds number based on diameter ranged from 
approximately 1.5x10° to 3x10°. This is on the edge of the transition region from laminar 
to turbulent separation, which is indicated by a rapid drop in the drag coefficient for a cir- 
cular cylinder in a cross-flow. Based on this, the flow about the model should have been 


transitional, and is shown to be so by the test results. The actual separation point during 
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Table 5.1 Numerical and Experiment Non-Dimensional Results 
the test should have been further aft than the calculated laminar separation case, leading to 
measured forces and moments somewhat less than the laminar predictions. That this is 


true lends some reassurance to the validity of the test results. 


Page 65 





Figures 28 and 29 show the laminar separation results along with the measured 
forces and moments for a, = 20°. The predicted normal force follows the measured trend, 
however, the predicted side force does not completely follow the experiment trend. The 
saine relationship holds for the moments: the predicted in-plane moment, M’, follows the 
experimental data trend, where the calculated out-of-plane moment, N’, has a very differ- 
ent characteristic. In both figures, the last data point seems to cause the problem. Only by 
“filling-in” the predicted curve, however, can meaningful conclusions about the results be 


made. 
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Figure 28. Predicted and Experiment Non-Dimensional Forces 
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Figure 29. Predicted and Experrment Non-Dimensional Moments 
5.2 Frequency Analysis 
The quality of the steady state force measured 1s related to the time varying effects 
present during the test. To get a handle on the vibration effects, a plot of the raw counts 


output for the normal force bridge vs rotation angle (or time) was created. The summation 


of 10 revolutions worth of data is shown in figure 30. The 20°, 30 ft/s run with the wing 
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attachment was considered the most probable for vibration effects. To avoid any convolu- 
tion of the result, the raw counts vice the processed forces were analyzed. As seen from 
the figure, there seems to be some scatter and possibly high frequency noise in the signal. 
The data points do follow the expected sinusoidal pattem, however, and when fit with a 


fourth order curve, produce a reasonable waveform. 
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Figure 30. Summed Raw Counts Output 
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A fast fourier transform (FFT) algonthm [17] was programmed and used to trans- 
form the 10 revolution raw counts output for both the N1 and N2 normal force channels. 
The FFT quickly accomplishes the discrete fourier transform of the N real data points 


(h,’s) as shown: 


The FFT produced a frequency and amplitude decomposition of the data points for fre- 


quencies up to the Nyquist or folding frequency. The Nyquist frequency is defined as: 
t= aA where A=sample time interval. 


For the 200 rpm runs, the folding frequency was 53.33 Hz. 

The FFT of the data used to produce figure 30 is shown in figure 31. The spike pres- 
ent at 3.3 Hz represents the fundamental rotation frequency (200 rpm). There are no major 
high frequency effects present, and only minimal lower frequency effects. This proves that 
the 3 Hz cutoff filters in the signal amplifiers did futer out any major high frequency noise, 
if any were present! For comparison, the inertia raw counts data was processed and is 
shown in figure 32. The figure does highlight the presence of some high and low fre- 
quency noise in the full-up data signal by comparison. The difference between the two 
figures could account for the scatter seen in the time varying signal. For further insight, 
the side force balance channels, Y1 and Y2, were also processed for the same full-up test 


run. The results shown in figure 33 are very simular to those for the normal force channels. 
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FFT of N1 & N2 Raw Counts 
200 RPM, 30 ft/s, 20 Deg 
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Figure 31. FFT of Normal Force Channels 
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FFT of N1 & N2 Row Counts 
200 RPM, in Air, 20 deg 
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Figure 32. FFT of Normal Force Channels, Inertia Run 
Because of the 3 Hz cutoff filters in the signal conditioners, no significant higher fre- 
quency effects could be seen (i.e., wall frequencies, higher harmonics, etc.). A look at the 
unfiltered shaft encoder output for rpm variation over a cycle however, could provide some 
high frequency information. Four cases were processed: 20°, 200 rpm inertia and full-up 
runs, and 20°, 25 rpm inertia and full-up runs. The cases covered the test spectrum for the 
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FFT of Y1 & Y2 Raw Counts 
200 RPM, 30 ft/s, 20 deg 
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Figure 33. FFT of Side Force Channels 
20° runs. The FFT results are shown in figures 34 and 35. The 200 rpm full-up nun defi- 
nitely shows a 2™ harmonic at 6.67 Hz. This is probably from the wing attachment passing 


by the support wing mounted to the forward nacelle. There are some other harmonics, 
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however, on the whole, the magnitudes are very small, typically tenths of rpm’s. No wall 
frequency component (4 x rotation frequency) was visible. Also, no 50 Hz component (the 


calculated torsional resonance frequency) was seen. 
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Figure 34. 200 RPM, Inertia and Full-up Run FFT Results 
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Figure 35. 25 RPM, Inertia and Full-up Run FFT Results 
RPM variation over a cycle is an umportant factor in obtaining accurate force and 
moment measurements. The FFT analysis produced the variation magnitude. The funda- 


mental rotation speed fluctuation was small for both rpm cases: 1.05 rpm for 200 rpm 
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(0.525%), and .18 rpm for 25 rpm (.72%) full-up runs. The 200 rpm output was much less 
than the calculated maximum rpm variation of 4.8 rpm. The low variance is a direct result 


of counter balancing the large sheave to account for the model/sting off-center weight. 
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Chapter 6 Conclusions and Recommendations 


6.1 Conclusions 


As part of continuing research into the flow about slender bodies of revolution, a 
coning motion apparatus for hydrodynamic model testing was built and demonstrated. The 
resulting test apparatus represents a culmination of effort by MIT, Stevens Institute of 
Technology, and General Dynamics, Electric Boat Division, and Convair Division. 

As the first test site for a coning motion experiment, the MIT variable pressure water 
tunnel required moderate modification. The modifications included installation of the 
model drive shaft, supports, and motor drive system. A workable data acquisition system, 
capable of capturing both the time varying and steady state effects, was established. The 


final system was capable of: 


. Positive and negative rotation rates to 200 RPM. 
° 0 to -20° coning angle settings in 2° increments. 
. Free stream velocities to 30 ft/s. 


The demonstration tests conducted covered a limited range of spin rates and coning 
angle settings. Reynolds number based on model length varied from 4.04x10° to 6.06x10°. 
Cross-flow Reynolds number based on body maximum diameter extended from 1.5x10° to 
3x10°, covering flow regimes in range of the transition from laminar to turbulent separa- 
tion. The steady state forces and moments were non-dimensionalized and presented as 


functions of the non-dimensional spin rate. Some of the main observations are: 


° Effects of Reynolds number variation are insignificant for the test range inves- 
tigated. 
° The out-of-plane forces and moments (Y’ and N’) are significant, reaching mag- 


nitudes of 30% of the in-plane force. 
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The forces and moments are non-linear functions of the spin rate for coning 
angle settings > 14°. 

The side force, Y’, is a possibly quadratic or higher order function of the coning 
angle. 

Measured test data correlates fairly well with numerical results, falling within 


the range of no separation to laminar separation. 


Overall, the first hydrodynamic captive model test using a rotary balance apparatus 


was a great success and opens the door to obtaining test data for motion in a non-planar 


cross flow. Ultimately, this will lead to a better understanding of the physics governing the 


behavior of freely moving bodies 1n large amplitude maneuvers. 


6.2 Recommendations 


Because of the extensive literature covering the experience of aeronautic researchers 


with the rotary balance apparatus, many first time mistakes were avoided. However, the 


different nature of testing in a water versus wind tunnel brought on its own host of prob- 


lems. The largest differences were the presence of buoyant forces, the thermal effects 


caused by immersion in water, and working in medium of roughly 1000 times greater 


density than air. For subsequent on testing, several recommendations can be made: 


Streamline the sector assembly to reduce both drag and cavitation caused by 
the sharp comers. 

Counter balance the apparatus at the model end to minimize rotational variance 
over acycle. Counter balancing 1s essential. 


Recalibrate the balance for a standard axis system convention. 
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Exchange the amplifier/filters for units with a higher cutoff frequency (if useful 
time varying information is desired). If this is not done, then include some 
separate monitoring system to check for vibration effects. 

For greater accuracy, conduct an on site calibration that accounts for local sys- 
tem peculiarities and determines model/sting deflection as a function of load- 
ing. 

Conduct all tests at a given @, in both rotation directions to account for flow 
angularity. 

Quantify the temperature sensitivity of the balance to allow for data correction 
later. 

Replace the 2 HP motor with a 5 HP motor. This will also minimize the RPM 


variance over acycle. The current motor is undersized. 
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Appendix A Taylor Series Expansion 
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APPENDIX A 


Taylor Series Expansion 


Al Taylor Series Expansion 


A Taylor series expansion with respect to the rectilinear and angular velocity compo- 
nents about a given condition is accomplished for representative forces and moments. The 
expansion will show the hydrodynamic derivatives that are of umportance and which relate 


directly to this study. 


A 1.1 Conventions 


The axis system, hydrodynamic forces, and the non-dimensionalizations are in accor- 
dance with reference [15]. The rectangular coordinate system in body fixed axis 1s as 
shown. The x-axis 1s chosen to coincide with the longitudinal axis of the body. The y-axis 
is in the honzontal plane 
and 1s positive to starboard. Body Fixed Axis System 
The x-z plane is taken to be 
the principle plane of sym- 


metry. 





Zw 
Z,N 
The hydrodynamic force and moments, X, Y, Z, and K, M, N respectively, and the 


linear and angular velocity components are non-dimensionalized as follows: 
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Taylor Series Expansion 


¢ Dimensionless Forces: 


x _  ¥ Ze 


pA? 1 AU? ~pAU? 











¢ Dimensionless Moments: 


mK /  M N 


ai 2 M =i 2 N 1 2 
1 pAIU 1 pAIU +pAlU 











« Dimensionless Velocities 


ae u Saal , WwW 
—— = 1 
U U U 
/ 
pl yall od 
U/ U U 
where: A = body frontal area 


/ = body length 
U = free stream velocity 


Note: 
2 Z 2 2 2 P 
U=u't+yvo 4+w lou" +v' +n" 
u’.v’,and wn" are direction cosines, not simply dimensionless u, v, and w, because 


the normalizing velocity, U, varies. 


A 1.2 Example Expansion of a Hydrodynamic Function 


Because of the tedious nature of the expansion, only the X hydrodynamic force is 
expanded as a representative force. The expansion is carried out about the condition of 


straight ahead motion at the instantaneous surge velocity, u(t): 
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Taylor Series Expansion 


u=u(t) vew=p=q=r=0 
or in non-dimensional form: 
fy res lah a sg) 34s) aac a0 
The X force is a function of the three linear and three angular velocities and can be written 


as: 
J 
Nee AI, Py gut) =5PA VaxRe aw. 4p <0 on ) Re =— 
Expanding about the reference point, the general form of the expansion is: 
aaa | 
riRev'on' prairie 5 All “2 + y’ Sp’ ged mer 2) 0, (Re,v, }, =, 7), C) 
where v, 0, €.7n. and € are dummy variables. 


Nee tien, 0.0,0,0,0) Re, =— 


Following reference [15], the standard shorthand notation is introduced: 





= -| 2x (Re,, ¥,0,0,0,0) 


v' sO 
Carrying out the expansion for the zero” and first order terms: 


, a ry? Gh ,O ae , oO win 4 
Ds al > tH ‘55 +P re ae pale (Re,,v,6,6,1, ¢) +... 


eee tN he tx Der x Gg + Xr | sae 
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Taylor Series Expansion 


Carrying out the expansion for second and some of the third order terms and employ- 


ing the shorthand notation: 
I 
X’ =X’, +[k = lterms] +7 [X/ e+ 2X" vw" + 2X4 vp! + 2X" vg! + OX vr’ 
+X" Ww + 2K" wp! + 2K" wg’ + 2X" wr’ 
, #2 , re , re 
+X op +2X pq +2X pr 
+X" .q” rs 2K" gr +X”) 
J >? re +? yo rf s2_ 9 , v2 , s2_ + 
+e lA Pe ero Neen be vaey eg + 3X00 5 


$30 we + GAG Pp OX in to. +X’ or D+... 
The expansion is the same for the other forces and moments: 
May 2 ak NG 
The expansion becomes quite cumbersome for a complete third order and higher 
expansion. Fortunately, many of the terms can be eliminated on the basis of symmetry. 


As a result of the x-z symmetry plane; X’,Z’,and M’ are even inv’, p’,r’ and 


Y".K’,and WN’ are oddinv’,p’,and_ r’. Using this and the data presented in Tables II - 


VI of Reference [10], the X and Y’ force terms are reduced to: 
] 
Ke tix +X g | eee a Ge a 
, 72 , Rear , #2 , roe , £2 , , 
Seed eget Pp Lk Dr +t Kg +X Fo) 


] 
ras = OCeT TS eee 
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Taylor Series Expansion 


ea y +[Y' vy + De +Y'r’ 


ee [Pree ry ov Git ly ip t 2Y ume + 2p g +2Y" qT | 


= ‘ 


] 
+—[k = 3terms] +.... 
6 
The force and moment terms can now be put in their final dimensional form. As a 


representative force, the X force becomes: 
re ] 2 re , J 
X=xX (3 pA us [Aw +X qf] pA v) 
] , 2 , ] 
+51% wY F2X  vwrit..] PA 


l ; 1{ pA 
ty [X’ vo 3X7 ew to] (2(%4)) +... 


The final form follows the same procedure for the other forces and moments (the dimen- 


sionalization ts slightly different for the moments K, M, and N). 


A 2 Data Analysis 


The Taylor Series expansion fits an assumed model of the form: 


Y’= Lof wi pig r)+e 
J 


where: Y° = measured hydrodynamic forces 
b, =coefficients 
jf, = fitting” functions 


€ = truncation error 
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Taylor Series Expansion 


The fitting functions are the controlled variables during the experiment (i.e., the result- 
ing linear and angular velocities and their combinations for a given rotation rate, coning 
angle, and tunnel velocity). The measured hydrodynamic forces and moments are the 
dependent variables resulting from the test runs. The basic matrix equation to solve is: 

{Y} = [X] {5} (1) 

where: { Y } = measured hydrodynamic force vector (1 x m) 

[X] = fitting functions (n x m) 

{b} = least squares fit coefficients (1 x n) 

m = # data runs 

n = # fitting functions 

{Y} and [X] are known frcm the experiment. The basic equation is unsolvable as 1s; 

therefore, a least squares approach 1s adopted. The method is to find the vector {b} that 


minimizes the sum of squares of the error. This is equivalent to minimizing the equation: 


([X] {b} — {¥})’ (x) {5} - {¥}) 
The minimization is found by multiplying both sides of equation (1) by [X]' and then 


solving for {b}: 
[X]" [X] {b} = [X]" {(¥} 


{b} =({X]" (XJ) 2X’ {¥} (2) 
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Taylor Series Expansion 


Once the coefficient vector {b} 1s known, it can be directly related to the corresponding 


hydrodynamic coefficients through a multiplicative constant. For example: 


X’=X',+bw'+b,q' +b,v" + 


b =X’. 
b aC 


2 hee 
2 wh 


5 
- 


b, = 


fe 
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Appendix B Inertia, Gravity and Buoyancy Forces 
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APPENDIX B 


B1 Inertial Forces 


The forces and moments experienced by the body-fixed balance due to the rotary 
motion of the apparatus are calculated as functions of the model moments of inertia (L,,, I,,, 
I,,), the rotation rate (@), and the coning angle (a@,). The approach will be completely gen- 
eral, and later adapted to this particular study. 

The axis convention 1s as follows: 

° Fixed (in Space) axis system: X,, Yo, Zo: 

. Body fixed axis system: x, y, Z. 


The coning angle, @,. , 1s defined as the angle between the x and x, axes. The x-axis is coli- 


near with the centerline of the model, the x,-axis is colinear with the centerline of the test 
section. The rotation rate vector, 0, is also colinear with the tunnel centerline. The 
nomenclature and axis system used here is consistent with the list of nomenclature. 

The equations of motion can be derived from a dynamic analysis. Using the basic 


Newtonian force/noment relationships: 


ey ee, 
Peg esteramtiatia)) Se ek Ue) (1) 
se 
Me= a (Angular Momentum) = 7 (Ho) (2) 


where : 


m = model mass 


Uc= velocity of model CG 


He = angular momentum of model at CG 
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For a reference point not at the CG, equations (1) and (2) can be rewritten as: 


a | 
F_=m—(U.+Q.xR 
c aha c c G) 
Bae «ee, ee 
M,=Me+Rex|“(mUe) 
dt 
where : 


Uc=uu,+vu,+wu,= Velocity of reference point in space. 


_——=3$ 


Q.= pu, +g u, +ru,= Angular velocity of body about the reference point. 
Rg =Xgl, t+ yg, +igt,= Vector distance of CG from reference point. 


Carrying out the differentiation and vector multiplication and setting the acceleration 


tenns = 0, the force and moment components reduce to: 
X =m[qu —rv ~x6(q° +1") + ¥o(Pq) +26 (Pr) 
Y =m[ru — pw — yor’ + p*)+2¢6(qr) + x6(4P)) 
Z = m|(pv — qu -z6(p° +9") +X¢(rP) + Yo(r9)] 
K =(/,—1,)qr +m[Ye(py — qu) — z¢g(ru — pw) 
+ X¢Ve(Pr) — Xeig(PG) + Yezg(r’- 9) 
M =(1,-1.)rp +m[zo(qw —Wv)-—x,6(pv -— qu) 
+ YoZg(QP) — YoXc(Gr) + Xezg(P - 1’) 
N=(1,-1,)pq +m[x_(ru -— pw) —ye(qw -rv) 
+ =gX¢("q) — zg Yo(P) + YorG(’ - PDI 


where X¢. Vg, 7g are the component distances of the model CG from the reference point. 
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For the formulation to be useful, the body-fixed velocities had to be defined with 
respect to the known variables, ® and a,.. A transformation matrix was derived that resolved 
the rotation rate (w), which was defined in the x,, y,, Z, coordinate system, into the body- 


fixed x, y,z system. The derived result (for no roll angle, = 0) is as shown: 


x cosa. sind®,sinof  —sina,cos@r]|| xX, 
Vist O Cos Wr Sin (wr y., 
4 sing. —cosa,sinox  cosa,cosm@r |} z, 


Given that the rotation vector is ®= @u,,, the relevant transformation is: 


p =(cos a.) 
Ae 
r =(sina@,)o 


Because the reference point is fixed, u =v =w=Q. Substituting into the force and moment 


equations: 
X =mw [-x,(sin’ o,) +z,(cos @, sin a,)] 
Y.=-m@ y, 
Z=m w’ [—z,(cos” Qa.) +x-(sin a, cos o,)] 
K =moy yo sina,[x, cosa, +z, sina, 
M =q@ [(l,-1) (sina, cosa.) +mXoz,(cos’ a, — sin’ «,)] 
N =-m@ yo cos @,[z, sin @, +X, cos O,] 
The 7 subscript refers to inertia forces. 


For the model used in the experiment, ye =z, =0 andi,= .40 inches. Because of the 


symmetry of the model about the x-z and x-y planes (both geometnic and inertial): 


° IL,=],=1, =0 


) 
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By inserting these relations into the inertial force and moment equations, the final form 


results: 


2 - 2 
X; =-MW xX, sin’ O, 


Y,=0 


2 : 
Z,=MW X_ SIN, cosa, 


M = (/,—/.) (sina, cos @,) 


K,=N,=0 





The forces and moments above are actually the reaction forces and moments. To be consis- 
tent with the inertial forces and moments measured by the balance, the negative of the 
relations should be used for comparison. Therefore, the negative of the relations were used 
to predict the inertia forces expenenced by the balance, and then compared with experiment 


results. 


B2 Gravitational and Buoyancy Forces 


Having derived the transfonnation matrix in the previous section, the calculation of 
gravitational and buoyancy forces (and moments) 1s straight forward. The gravity and buoy- 


ancy vectors are: G= Wu, and B=-B u,,. The forces follow directly: 
Xy_p =—-(W —-B)sina, COS GW 
Yy_, =(W —B)sin or 


Zy 3 =(W -B)cosa, cos wt 
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For the moments, a simple moment balance at the reference point results in the relations for 
weight and buoyancy: 


K,.=W[v, cos @, cos Wf — Zz, sin a] 
M,.=—-W cos of[z,.sin®, +x, cos & J 


Ny = W[%, sin of + yo sin &, cos Wr] 


The same equations result for buoyancy by substituting -B for W and Xp, yg, 2, for XG, Ve. ZG. 
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Appendix C Test Data 


The following pages contain the steady state inertia and test run data for the plots pres- 


ented in Chapters 3 and 4. 


Page 96 


. al 
nt her oS AD 


J ese? eee 


iG 


~ ne ’ 4 





C1 Inertia Data 


qa = -20° 


Cc 


Measured Forces & Moments SNAME Axis Convention 
Calculated Forces & Moments 
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Measured Forces & Moments SNAME Axis Convention 
Calculated Forces & Moments 
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C2 Water Test Data 
= .20°, V = 20 ft/s 


Dimensional Forces and Moments: 
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a, = -20°, V = 20 ft/s 


Non-Dimensional Force and Moments (SNAME Axis Convention) 
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oO, = -20°, V = 28 ft/s 


Dimensional Forces and Moments: 


po tv{| x | y fz] xe | ™ | on 
a ee 
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a, = -20°, V = 25 ft/s 


Non-Dimensional Force and Moments (SNAME Axis Convention) 
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a, = -20, V = 30 ft/s 


Dimensional Forces and Moments: 
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a, = -20°, V = 30 ft/s 


Non-Dimensional Force and Moments (SNAME Axis Convention) 


~ 9.0000! 0.12469! 0.039411 0.90976! 0.00013] -0.27845! -0.01913 
0.1709 0.13105! 0.96614] 0.00024! -0.28222] -0.05179 
0.20240 -0.00008} -0.26908| -0.09011 
0.129131  0.25327/ 1.057131 0.00007] -0.25428] -0.10780 
0.28506 0.00002) -0.24636| -0.12362 
0.6888 0.122171 0.35217] 1.19665} -0.00010] -0.22299| -0.14844 
0.7735 0.38372 0.000081 -0.21574| -0.15882 
0.8595| 0.12194 0.000111 -0.20876! -0.16937 
0.11778 0.00013! -0.20177| -0.17490 
1.2200! 0.09775 1.56073! 0.00015] -0.17097] -0.19499 
0.09364] 0.46502 0.00039! -0.16090| -0.19673 
0.08902 0.00061] -0.15091] -0.19670 
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a. =-14', V = 25 ft/s, Positive Rotation 


Dimensional Forces and Moments: 


ete ff 2 
A 
Tea] asco] ass] ose] tiara] oanis| viesre] —saiay 
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a. =-14°, V = 25 ft/s, Positive Rotation 


Non-Dimensional Force and Moments (SNAME Axis Convention) 


Pe To Tr Te Te Tw To 
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a. =-14°, V = 25 ft/s, Negative Rotation 


Dimensional Forces and Moments: 


ef tt 
coal asa] 2am] onal Tose aan] vost] asi 
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a. =-14°, V = 25 ft/s, Negative Rotation 


Non-Dimensional Force and Moments (SNAME Axis Convention) 


i 
T08Ts 
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i = -8 , V = 25ift/s 


Dimensional Forces and Moments: 


2 i} Se SS 
Toca] 2500] —19aea] -a07e7] —a40os] -oa7e0] —aoxala] sean 
Tiresilasoi] —aa0es] areas) cario] -oas0o) eaooee] waste 
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oa, = -8°, V = 25 ft/s 


Non-Dimensional Force and Moments (SNAME Axis Convention) 



















aoooer| ones] ones 
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Appendix D FORTRAN Routines 


The following pages contain the FORTRAN routines used for obtaining the data pres- 


ented in the study. Also included is the FFT algorithm used for the data analysis. 
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$DEBUG 
program DC 


OO AO YOY a ae ae 


a mee 


OO os oo) OO Ore OY OO EY Or era | 


‘o) 


Wea OO Oa ) Oe 


Version: 1.2 
Date: 03 April 1989 


Purpose: Main driver program that links with the following 
programs: Davlab, Common, DCDave and Calibr.For. 
This program 1s to obtain readings from the balance, 
store them (ac and dc), and convert the steady state 
components to forces and moments. 
Program for ROTATING measurements, DC forces 
and moments. 


Programmer: Dave Johnson and Glenn Mckee 
Language: Microsoft Fortran 4.1 


Variables: 
fdata(i): Raw MIT counts array filled with counts 
from al] a/d channels. ("full” data) 


pdata(i): Raw MIT counts array with 6 preferred 
channels of a/d channels. Array values 
are later corrected for zeros and convair 
/MIT counts. This array 1s then mult. 
by a(i,j) for the balance loads. 
(‘preferred” data) 


loads(i): Array with balance loads. loads(i) comes 
from the matrix multiplication of pdata(i) 
and a(i,j). 


zerot(i): Row array with averaged counts for all 
a/d channels. This array is read by the 
subroutine "filter" to convert to the 
preferred 6 channels and put into zeros(1). 


a(i,j): Coefficient matrix [B]. From {L} = [B]{R} 
Files: 
mame: Raw MIT counts file with data from all 
a/d channels 


Ratio: file of Convair R-cal/MIT R-cal readings for the 
6 preferred channels. 


M2X1T.DAT: Balance coefficients for Roll bridge 2 and 
Axial bridge 1, referenced to C of Rot. 


*** Data Conventions Used in This Code *** 
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OOS OO AO ae ae Fe 


Index Readings() Loads() 
] Normal Force 1} Normal Force at center 
Z Normal Force 2 Pitch Moment at center 
3 Side Force 1} Side Force at center 
4 Side Force 2 Yaw Moment at center 
> Selected Roll Bridge Roll Moment at center 
6 Selected Axial Bridge Axial Force at center 


( The zeros() and ratios() vectors use the same conventions as 
the pdata() ) 


real pdata(6), loads(6), zeros(6), ratios(6), fdata(9) 
real zerot(9), dpc, rpmc, dp, rpma,iloads(6) 
character*3 tcase, zcase 

character*!] zname, Iname, mame.iname.calibr 
character* 8&0 id, header,head! .headc 

character ans 


common /dpcell/ rho, zero, full, tt 
common /rpm/ rvz, rvcal 
common /revsect/ nrev.nsect 


-- open files 


open (unit=10, file = ‘ratio’ , status = ’old’ ) 


write (*,*) 
smite (*,") ° *** PROGRAM DC.FOR *** °’ 
write (*,*) 


cdb --- Debug commands for reading in Convair test case. This 1s to 


cdb 
cdb 


check the data reduction algorithm ...... 


goto 52 
write (*,*) ' Note: Zeros for this program are stored in the’ 
emite(™,*) file zname. The program zero should be’ 
evmite(*,*) run first to get the proper DAILY zeros. ’ 
write (*,*) 


write (*, °(A\)’ )’ Input the zero (counts) file name: ’ 
read (*,1) zname 
write (*,*) 


-- Read in zeros from file zname 


open (unit = 2, file = zname , status = ’old’) 

read (2,’ (a80)’) header 

read (2,’ (a80)’) id 

read (unit =2, fmt = 2, err = 51) zcase, (zerot(1) , i = 1,9) 
write (*,*) 

write (*,*) ’ Zeroes read in successfully! ’ 

write (*,*) 


--Debug command 


write (*,*) zerot(] ).zerot(9) 





goto 52 


oH write (*,*) 
write (*,*) ’ Error in reading zero file ’ 
goto 100 

DZ continue 


C ----- Read in ratios from file of preferred channels (file "RATIO") 
read (unit = 10, fmt = 6, err = 5S) (ratios(i) , 1=1,6) 
write (*,*) 
write (*,*) ’ Ratios read in successfully’ 


C ----- Debug command 
write (*,*) ratios(1),ratios(6) 
goto 56 

i.) write (*,*) ° Error in reading in ratios’ 
goto 100 

56 continue 


cdb --- Debug for reading in Convair test case .... 
cdb goto 57 


C ----- Call filter subroutine for producing zeros(1) 
call filtz (zerot,zeros) 


ai/ continue 

C ... assign a vacant FORTRAN unit number ... 
lunit = 3 

c .. use the CONVAIR preferred set of coefficients (R2X1) ... 
item = 4 


call intcoef( Junit. item ) 
Write (*,*) ’ Coefficient Data has been read ° 


C ... these are the A/D readings without applied extemal forces ... 
call setzero( zeros ) 


5 ... these are the ratios of the current signals divided by the ... 
C ... Original signal levels at calibration ... 
call setrcal( ratios ) 


cdb --- Debug for reading in Convazrr test case ...... 


cdb goto 58 

C ---- Offer option of either taking data or processing a previous 
C run's data..... 

70 write (*,*) 


write (*,*) ’ Do you wish to take new data or process’ 
write (*,’ (a\)’) © old data ? (N/O): ’ 

read (*, fmt=9) ans 

if (ans .eq. ’n’) then 


C ----- Calibrate or read from an old calibration file? 
write (*,*) 
write (*,’ (a\)’) ’ Calibrate? (y/n):’ 
read (*,9) ans 


ae 





if (ans .eq. ‘n’) then 
write (*,*) 
write (*,’ (a\)’) ’ Input old calibration file:’ 
read (*,1) calibr 
open(unit=5, file=calibr,status=’old’) 
read (5,’ (a80)’) headc 
read (5,’ (a80)’) id 
read (5,15) rho,zero,full,tt,rvz,rvcal 
close (unit=5) 
goto 80 

endif 


C ----- Call calibration routines ..... 
call dpcal 
call rpmcal 


C ----- write result of recent calibration to file Calibr.... 
write (*,*) 
write (*,’ (a\)’) ’ Input new calibration file:’ 
read (*,1) calibr 
write(*,*) ° Type in new calibration file header:’ 
read (*,’ (a80)’) headc 
open (unit=5 ,file=calibr) 
write (5, (a80)’) headc 
write (5,17) 
write (S,15) rho,zero, full,tt.rvz.rveal 
close(unit=5 ) 


80 write (* ,*) 

o---- write (*,’ (a\)’) ’ Input # of revolutions, nrev:’ 
C ----- read (*,3) nrev 

@e----- Merite (4) 


write (*,’ (a\)’) ’ Type in tcase (3 max): ’ 
read (*,’ (a3)') tcase 


miseet = 32 
mey — 10 
C nreeeee-- Go to data taking subroutine 


Call takedata (fdata.dpc.rpmc,tcase,header) 
write (*,*) ’ takedata finished! ' 


rite (* 
C --------- debug command 
cdb pause 

goto 62 

endif 

cdb --- Debug for reading in Convair test case or for inputing 
c previous run file .... 
58 continue 


mrite (*,*) 

write (*,’ (A\)’) ’ Input PREVIOUS (DC) raw counts file: ’ 
read (*,1) mame 

write (*,*) 


Blise 





open (unit=9, file = mame, status = ’old’) 

read(9,’ (a80)’) header 

read(9,’ (a80)’) id 

read (unit=9,fmt=1 8 ,err=61) tcase,(fdata(i), i=1,9),dpc,rpmc 
write (*,*) ’ Previous counts ..... 
write (*,*) fdata(1), mmc 


C ----- Read in old calibration file............ 
write (*,*) 
write (*,’ (a\)’) ’ Input old calibration file:’ 
read (*,1) calibr 
open(unit=5, file=calibr,status=’old’) 
read (5,’ (a8Q)’) header 
read (5,’ (a80)’) id 
read (5,15) rho,zero,full,tt,rvz,rvcal 
close (unit=5 ) 


cdb do 60j=1,6 


cdb zeros(})=0.0 

cdb60 continue 
goto 62 

61 write (*,*) ' Error in reading Previous counts....’ 
goto 100 

62 continue 


‘es 3H she ae he he Cae ae IC ae ae ie ke ake ake ak ae ae 2k ake oc Data Reduction Portion aie he ak ae hae ake ake RC akc ake ake ae ake ak ak ake ake akc ake ake akc ak afc 
C ----- Call filter subroutine to reduce the readings to the preferred 
c 6 channels (1.e. read 1n data to pdata(1)) 

call fulter(fdata, pdata) 


write (*,*) ° filter done ’ 


C ----- debug command 
cdb pause 
merite (* ,*) 


write (*,*) ° Now Processing Data!!! ’ 
write (*,*) 


C ... adjust the values by removing the zeros and compensating .... 
C ... for the ratio of the signal level ... 
call doadjus( pdata ) 


call doload( pdata, loads, ier ) 
if( ierr .ne. O ) then 
write(6,12) tcase 
endif 
20 continue 


C ----- Convert dpc and rpmc counts to velocity and rpm... 
call rpmvel(dpc.rpmc,dp.rpma) 


a. Write loads to screen 
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write (*,*) 

write (*,*) ’ Loads for this mun... ’ 

write (*,*) 

write (*,10) 

write (*,7) tcase, (loads(1) , 1=1,6),dp,rpma 


¢ ---- Ask if INERTIA run or HYDRO run. This determines whether or not 
c _ the inertial forces are subtracted out from the above loads..... 
75 write (*,*) 
write (*,’ (a\)’) ’ Is this an inertial or hydro run? (i/h):’ 
read (*,9) ans 
if (ans .eq. *h’) then 
write (*,*) 
write (*,’ (a\)’) ’ Input Inertial FORCE & MOMENT file:’ 
read (*,” (all)’) iname 
open (unit=6, err=75 ,file=iname,status=’old’) 
read(6,’ (a80)’) head] 
read(6,’ (a80)') id 
read (6.16) tcase,(loads(1), 1=1,6) 
close(unit=6) 


C --------- Do Load subtraction ........ 
qos) i=) ,6 
loads(i )=loads(i)-Uloads(1) 
&5 continue 
endif 
C ----- Wnite loads to file Iname......... 


write (*, ' (A\) ° )’ Input the load output file name: ’ 
read (*,1) Iname 
write (*,*) 


open(unit=7, file = Iname, status = ’new’ ) 
C ----- write (*,*) ’ Input header for load file... .’ 
C ----- read (*,’ (a80)’) header 

write (7,’ (a80)°) header 

write (7,10) 

write(7,/) tcase, (loads(1), 1=1,6),dp,rpma 


c ---- Write loads to screen again.......... 
write (*,*) 
write (*,*) ’ Final loads for this run... ’ 
write (*,*) 
write (*,10) 
write (*,7) tcase, (loads(i) , i=1,6),dp,mpma 


write (*,*) 
90) write (*,’ (A\)’” )’ Would you like another test run? (y/n) : ’ 
read (*,fmt =9, err = 90) ans 
m=m+1 
if (ans .EQ. ’y’) goto 70 


C ----- omnat Statements .....-......c:.00--..0+ 
] format(Al11) 


alle = 





Z format (A3, 9F7.1) 

3 format (12) 

6 format (6F8.4) 

7 format (A3, 6F9.4,2f9.2) 

9 format (A1) 

10 format(’TC’ ,4X,’Z’ ,8.X,°M’ 8X,’ Y’,8X,’N’,8X,’K’, 
+ 8X, X’,8x,’ VEL’ ,6x,"RPM’) 

12 format(’ Convergence problem on point ’,A10) 

1S format (6f9.5) 

16 format(a3 ,6f9.4) 

17 format (3x,’Rho’,6x,’Zero’,5x,’Full’,6x,’Tt’,6x,’Rvz’, 
- 6x,’Rvcal’) 

18 format (a3,11F7.1) 


C ----- Close all remaining open files ..... 
close(unit = 10) 
close(unit = 7) 
close(unit = 1) 
close(unit = 2) 


100 continue 
end 
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$DEBUG 


C anna nn en ee ne-e------------ PROGRAM DCDAVE.FOR --------------------- 
c 

C ----- Version: 1.2 

C 

C ----- Note: Full Revolution Version , Therefore nrev=10 max! 

c 

C ----- Date: 29 March 1989 

C 

C ----- Purpose : To provide the collection of subroutines necessary 
c for rotating data taking and DC processing. This 

Cc program is for linking with the main driver program 

C DC.FOR. 

C 

C ----- Programmer: Dave Johnson 

c 

C ----- Language: Microsoft Fortran 4.] 

C 

C ----- Subroutines: 

C TAKEDATA: Primary rotary data collection sub- 
C routine 

C 

C FILTER: Reduces readings to the preferred 6 

C channels for later processing 

C 

c Rien Reduces zero readings to 6 preferred 

C chanels for later processing 

C 

c RPMVEL: Converts dp and rpm counts to actual 

G tunnel velocity and model rpm 

e 2k ke ae She ake ake ote ie ake ote Sie ake ote ae ae aie ae af ae ake af abe fe afc ate ake af ape ake fe ake ate abe ak ake ate ae ake oe ae aie ake ae oe ake ate aie ake oe aie ate ake aie af ake aie ade af oe af ae ate af ate ae ae ae ake ote 


SUBROUTINE TAKEDATA(fdata,dpc ,mpmc,tcase header) 


C weene-n------- 22 2--2------- SUBROUTINE TAKEDATA---------------------- 
C 

C ----- Arrays: _ fdata(i,j): Averaged (over nsamp) raw MIT counts 
C for all 9 a/d channels of balance data. 

S 

c larray(i): Row array filled by EACH mcbatod call 
C 

C tarray(i,)): Temporary array with rows of iarray 

é and columns of nrev*nsect samples 

e 

€ bigarray(i,j): Summed counts array that 1s divided 

G by nsamp and written to tares(1,}) 

C 

¢ dpc(i): Raw MIT counts read from dpcell 

C 

c rpmc(i): Raw MIT counts read from daytronix chan C 
C 

C ----- Files: RCxxxx.xxx: Raw data file in MIT counts. 

c 


- 119- 





----- Dimensions set for max of 32 sectors/rev x 10 revs = 96 data pts 

and max of 1] channels sampled at 3 samples/channel ...... 

integer* 2 iarray(33), tarray(320,33), nval, n, m, nsamp 

integer*2 nchan, echan, k 

real rfdata(320,9), bigarray(320,11), frpmc(320), fdpc(320) 

real rpmc,dpc,rpma,dp,sum 1] (32,11),fdata(9) 

character* 11 mame,rsum,tsum 

character*80 header 

character*3 tcase 


common /revsect/ nrev, nsect 


nenan = || 
echan = 10 


----- Setclock set for 1kHz 
call mbopen 
call setclock(.001,1) 


----- Input data counts filename and # of samples to be taken at 
each data point ..... 

mite (**) 
write (*,*) ’ Format for Raw Cts File is RCXxXxxx.xxx ’ 
write (*,*) ’ where xxxxx.xxx Is tst case, V, alpha, rpm ’ 
smite (*.* } 
write (*, ’ (A\)’)’ Input Raw Counts Output Filename : 
read (*, 5) mame 
open (unit = 1, fue = mame) 


9 


write (*.*) 
----- write (*.*) ’ Before entering # of samples per data point, ’ 
wo--- write (*,*) ’ remember that for any high speed data taking, ’ 
= write (*.*) ’ nsamp should not be > 1. ' 
aoe write (*,*) 
~---- write (*, ’ (A\) ’)’ Input # of samples for each data pount: ’ 
o---- read (*, 6) nsamp 

nsamp = | 


write (* ,*) 
pause ’Press any key to start data taking ....... 
----- The following logic is based on inputs to the digital input 
ports #0 and #1 being 32(J) and 1/rev(K) pulses respectively. 


The elaborate if-then statements ensure that the data taking 
actually starts at 0 degrees and counts for nrev revolutions. 


——= Initialize counters ...... 


nm 0 
m= 0 
k=0 


call dina(nval) 
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C ----- If the following expression is true, then K is low.... 
100 if (nval .le. 1) goto 200 
cal] dina(nval) 
goto 100 
C ----- If the following expression 1s true, then K is high and thus 
C the model is at 0 degrees. Now start counting 32 times for 1 
c revolution and then start taking data on 33rd count...... 
200 call dina(nval) 
250 if (nval .ge. 2) goto 290 
call dina(nval) 
goto 250 
C ----- Start counting 32 times off J trigger 
290 call dina(nval) 
300 if ((mod(nval,2)) .EQ. 1) goto 400 
call dina(nval) 
goto 300 
400 n=n+]1 
if (n .EQ. 32) goto 500 
C ----- Look for J low ..... 
450 if ((mod(nval,2)) .EQ. 0) goto 290 
call dina(nval) 
goto 450 
C ----- Look for J low before taking data .... 
500 if ((mod(nval.2)) .EQ. 0) goto 600 
cal] dina(nval) 
goto 500 
C ----- Look for J high to trigger data taking. This ensures start at 
c OTGeSrees.... 
600 call dina(nval) 
650 if ((mod(nval,2)) .EQ. 1) goto 700 
call dina(nval) 
goto 650 
C ----- Now at 0 degrees, begin taking data ........ 
700 call mcbatod (0,echan,1,nsamp, larray) 
m= m+1 
C ----- Throw temporary data into tarray to prevent overwrite of data .... 
do 750 i=1, (nsamp*nchan) 
tarray (m,1) = larray (1) 
750 continue 
C ----- Exit condition; exit loop after nrev revs of data taking (10 max) 
if (m .EQ. (nrev*nsect)) goto 755 
C ----- Look for J low, if J is low, go up to top of loop and look for 
C J high again. This prevents multiple tnggers off the same 
e pulse, especially important for slow rpm runs ....... 
Fol if ((mod(nval,2)) .EQ. 0) goto 600 
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82 


call dina(nval) 
goto 75] 


-- Done taking data, now put into proper arrays and store in file... 


write (*,*) 
write (*,*) ’ Done taking data ..... 
write (*,*) 


-- Fill array bigarray(1,j) with zeroes .... 


do 760 i= 1, (nrev*nsect) 
do 760 j = 1, (nchan) 
bigarray(i,j) = 0.0 
continue 


-- Place sum of nsamp number of data points in each storage bin .... 


do 780 1 = 1, (nrev*nsect) 
do 780 j = 1, (nchan) 
bigarray(1,}) = bigarray(1,}) + tarray(i, G+(nchan*k))) 
continue 


k = k+]1 
if (k .EQ. nsamp) goto 78] 
goto 765 

continue 


- Next, divide the samples by nsamp and put in arrays rfdata(1,j), 
fdpc(i) and frpmc() 
do 791 1= 1, (nrev*nsect) 
do 790j;=1,9 
rfdata(i,}) = bigarray(1,})/nsamp 
continue 
fdpc(i) = (bigarray(i,10))/nsamp 
frpmc(1) = (bigarray(1,11))/nsamp 
continue 


-- Finally. put data values into the data file mame ...... 


write (*,* ) 
mete (",~) ype in file Header ....... 
read (*,’ (a80)’) header 
write (1,’ (a80)’) header 
write (1,14) 
do 795 1= 1, (nrev*nsect) 
write (1,15) (rfdata(i,j) , j = 1,9), fdpc(i), frpme(i) 
continue 


-- Sum up nrev points into | rev of data and put into file ... 


Fill sum1, fdata, dpc, and rpmc with zeros .... 


goo2i=i,11 
do 82 j= 1, nsect 
sum] (7,1) = 0.0 
continue 
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90 


do 811=1,9 
fdata(i)=0.0 
continue 

dpc = 0.0 

rpmc = 0.0 


-- First do balance counts .... 


k=0 
do 85i=1,9 
do 85) = 1,nsect 
sum 1 (j,i) = sum1G,1) + rfdata(j+(nsect*k),i) 
continue 
k=k+1 
if (k .eq. nrev) goto 86 
goto 83 
continue 


-- Now do dp and rpm measurements ..... 


k=0 


continue 
do 90 j = I,nsect 
sum1(j,10) = sum1Q.10) +fdpcg+(nsect*k)) 
sum] (j.11) = sum1Qj,11) + frpmcGg+(nsect*k)) 
continue 
k=k+1] 
ik eq. nrev) goto 9] 
goto 89 


-- Divide by nrev for averaged values over 1 revolution ... 


do 1101=1, 11 
do 110) = 1, nsect 
sum1(j,1) = sum1(j,1)/nrev 
continue 


-- Now sum up over remaining revolution for averaged counts ..... 


do 951=1,9 
do 95 j= 1, nsect 
fdata(i) = fdata(1) +sum] Gj,1) 
continue 


-- Again, do dp and rpm .... 


do 105 j= 1 , nsect 
dpc = dpc + sum14,10) 
rpmc = rpmc + sum1Q@j,11) 
continue 


-- Divide by nsect ...... 


do 120 1=1,9 
fdata(1) = fdata(i)/nsect 
continue 


dpc = dpc/nsect 
rpmc=rpmc/nsect 


Re 





C ----- Wnite sum] and fdata, dpc, rpmc, to files ... 
write (*,*) 
write (*, ’ (A\)’) ’ Input summed (over | rev) file: ’ 
read(*,1) rsum 
write (*,*) 
write (*,” (A\)’) ’ Input TOTAL summed counts output file: ’ 
read (*,1) tsum 
open (unit = 15, file = rsum, status = ‘new’) 
open (unit = 16, file = tsum, status = ’new’) 
write (15,’ (a80)’) header 
write (16,’ (a80)’) header 
write (15,14) 
write (16,10) 
write (16,7) tcase,(fdata(i), 1 = 1,9),dpc,ppmec 
do 130 j= 1,nsect 
write(15,15) (sum1 (j,i), 1 = 1,11) 
130 continue 


C ----- FORMAT statements ............. 
] format (al 1) 

5 format (al 1) 

6 format (12) 

7 format (a3,11f7.1) 

10 femman 1C 4x NI’ 5x, N2’ 5x, Y1’,5x, Y2',5x, RI’, 


+ 5x,’ R2’,5x,’X1",5x,’X2’,5x,’X3°,5x,’DP’ ,4x,’RPM’) 
14 Merman sax. NI 5x. N2 5x, Y1 5x, Y2 5x, R1, 
+ Oe OX ol OX, ek AS BOK, DP 5x, RPM’) 
1S format (11F7.1) 
1000 continue 
close (unit = 1) 
retum 
end 


c (222 ff F See Se SSS ES SSS ESS SESS ESE SES SS SESE SESE SSE SESE SES EES SES ES SE SE 


subroutine filtz(zerot.zeros) 


C ------------------------ SUBROUTINE FILTZ --------------------- 
real zerot(9), zeros(6) 
ao 50 i= 1,4 
zeros(1) = zerot(1) 
30 continue 


zeros(5) = zerot(6) 
zeros(6) = zerot(7) 


retum 
end 


CFIC II OGIO GIGI IOI GI GI IOI GIGI II III IOI Kak iC ak aI ak IOI ICICI kK 


subroutine filter(fdata, pdata) 
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c 
C ----- Date: 16 Feb 1989 
C 
C ----- Purpose: Filter subroutine to take all channels of data and read 
c in the 6 preferred channels for processing by Maintst. 
€ 
C ----- Vaniables : 
c 
c fdata(i): array with all channels of data (11) 
c pdata(i): array with 6 preferred channels of 
C balance data. 
c 
C wee w wn www ww ww en on nn en wn ce ene ce cee msec ences sce ec cces esses cesssrrre= 

real pdata(6),fdata(9) 
C ----- debug command 
cdb write (*,*) ’ Have reached subroutine filter ’ 
cdb pause 
C ----- Do test reading conversion 

do 101=1,4 

pdata(i) = fdata(1) 

10 continue 
C ----- debug command 
cdh write (*,*) ° First 4 points read in successfully ’ 


pdata(5) = fdata(6) 
pdata(6)= fdata(7) 


retum 
end 


CO IOI Oli ii Iii i IOIOI IG IOI GIO iii IOI IGRI I II GIOIOIGIGI I IGIOIGG III III GIGIGI III I GIGI GIGI OIG IO iO Ok 


subroutine rpmvel (dpc, rpmc, dp, rpma) 


C --------------------------- SUBROUTINE RPMVEL --------------------- 
C 

C ----- Purpose: To convert dpcell and rpm readings to actual 

C velocities and rpm’s. 

c 

CC a re ee rn wenn ee eee eee ee ee 


real dpc, rpmc, dp, rpma 
common /dpcell/ rho,zero,full tt 
common /rpm/ rvz,rvcal 


C ----- Convert dp counts to voltages and convert to ft/sec 
fe using 409.6 counts/volt and Lisa’s function dpspeed ........ 


dpc = (dpc)/409.6 
dp = dpspeed(dpc) 


- 125 - 





----- Convert rpm counts to rpm's 
rpmc = (rpmc)/409.6 
Ipma = (revs(rpmc))*60/512 


retum 
end 
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real pdata(6), loads(6) 
integer ierr 


*** The subroutine takes the readings corrected for zeros and *** 
and scaling, and then converts them to equivalent loads 
using a linear 6th order fit determined by the Davidson 
Laboratory MCB program. 


She fe ae Be ae ake ae ake ae ake ake ake ae ae ake ae ake ake ape ae ae ale ae he ae akc ae he ae ae abe abe ae Bc aie ale Be She ae She she he ae he He she She he eae Be he he he he ae ahs ae ah ae ae ae ae ae ae 


common /convatir/ a(6,6) 


... there 1s no possible error in this case ... 


ierr = 0 


... do the matrix multiplication ... 


do 20 1=1,6 


loads(i) = 0.0 
do 10j;=1,6 
loads(i) = loads(1) + aGU,1)*pdata(y) 
continue 
continue 


... done, so returm ... 


retum 
end 


subroutine intcoef( runit, item ) 
integer runit, item 


*** Subroutine to obtain the necessary information from a file *** 
and to store it for later use of the conversion routines. 


RUNIT: An available FORTRAN unit number for reading the 


data file 
ITEM: The set of coefficients selected for use in the 
conversion routine. 
ITEM # Selection 
Roll Bridge Axial Bridge 
] l l 
2 l 2 
) ] 3 
4 p l 
5 2 2 
6 2 3 


BK ake ake ok ak ak ak ai a ak ak IC ak akc ake ake ak ake ake ake ake ake ake ake ake ake ake ake ake ake ake ak ake ake ake ake ake ake ak ake ake ake ak ae ake ak ake ake ake oe ak ake ake ake ake ake ake ake ake ak ake ake ake ake ae of 
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character*80 id 
character*9 name(6) 
common /convair/ a(6,6) 


data name/ ’M1IX1T.DAT"’, ’M1X2T.DAT’, ’M1X3T.DAT"’, 
"M2X1T.DAT’, "M2X2T.DAT’, ’M2X3T.DAT’/ 


*** attempt to open the data file *** 
open( unit=runit, file = name(item), status = ‘old’, err= 70 ) 


*** insure that zero is the default value *** 
do 10i=1,6 
do 10 j= 1,6 
a(i,j) = 0.0 
continue 


*** read the contents of the file *** 


read(runit.15) id 


format(a80) 


do 20i1=1,6 


read(runit,* end=25,err=25) ( a(i,j), j=1, 6 ) 
continue 


goto 30 


print *,, ERROR - Trouble reading Coefficient Data’ 


retum 


continue 
write(6,40) id 


debug commands 
open (unit = 6, file = ‘loadcoef’, status = ’new’ ) 
format(1x,a80) 
do 60:1=1,6 
write(6,50) 1, (a(i,j), j=1,6) 
format(1x,14,6f12.6) 
continue 
debug command 
close (unit = 6) 


close( runit ) 
retum 


wnite(6.80) name(item) 
format(5x,’ ERROR - Data file ’,a9,’ can not be read’/ 
aK. All coefficients are zero’//) 


end 
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won n--------------------- Subroutine doadjus------------------------------ 
subroutine doadjus( pdata ) 
real pdata(6) 


*** This subroutine takes the current readings, corrects for the *** 
readings for zero applied load, and then scales the result 
back to the levels at which the calibration was done. The 
retumed values should be the bits that would be read at 
CONVAIR calibration bench. 


(Peete SESS SES SEES ES EEE ESE EEE SEES ES ESSE ESSE EE EEE EEE EE EE EE ESS SE SE SE SS 


common /convar1/ zero(6), rcal(6) 


do 101=1,6 
... remove the zero first ... 
pdata(1) = pdata(1) - zero(i) 
pdata(i) = pdata(1)*rcal(1) 
10 continue 


retum 
end 


subroutine setzero( zeros ) 
real zeros(6) 


*** This subroutine takes zero applied load readings and stores *** 
== them for later use *** 


common /convar!/ zero(6), rcal(6) 
do 10i1=1,6 
zero(i) = zeros(1) 


10 continue 


retum 
end 


subroutine setrcal( ratios ) 
real ratios(6) 


*** This subroutine takes calibration signal ratios and stores *** 
=. them for later use *** 


common /convar1/ zero(6), rcal(6) 
do 101=1,6 


rcal(1) = ratios(1) 
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$DEBUG 
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10 


100 
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- This program contains the calibration subroutines for calibrating 
the dpcell and rpm channel. 


SUBROUTINE DPCAL 


Calibrates the differential pressure cell. 
Stores the calibrations "zero", and "full" 
in the common block called “dpcell". 


After calling DPCAL once in your program, 
a special function exists in the MHL library 
that you can use to convert voltages from 
the Daytronics channel C to velocity. It 

is called “dpspeed”, and has one calling 
parameter, ‘voltage’. 


Example: 


call dpcal 


<do an A/D conversion on channel 2 (Daytronics C)> 
freestream=dpspeed(voltage) 


> ete 


You do not have to include the common block in your 
main program to make these routines work 


Lisa Shields 
June 12, 1987 


common /dpcell/ rho,zero, full, tt 
integer*2 channel,ig.ival 
character ans, line*70 


write(* ,*) 
write(*,*) ’ Welcome to Dpcal!’ 
write(* ,*) 
write(*,’(A\)’) ’ Please input water temp (deg f): ’ 
read(*,* ,err=10) tt 
rho=1.9574-0.00028*tt 


FORMAT(5x,a70) 
write(*,*) 
write(*,*) 
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line=’The dpcell should be connected to the Daytronic channel c’ 
write(*,100) line 
write(*,*) 
line=’The span and balance should be set such that the digital’ 
write(*,100) line 
line=’display reads from -820 (zero flow) to +600 (with ° 
write(*,100) line 
line=’the "-cal" button depressed.) ’ 
write(*.100) line 
write(*,*) 
write(*,*) 
write(*,*) 
write(*,101) 
101 format(5x,’Please hit return after checking this ’,\) 
read(*,104) ans 
write(* ,*) 
line=’Sampling zero flow voltage for 5 seconds and averaging’ 
write(*,100) line 


call mbopen 

call setclock(0.01,0) 

zero=0.0 

channel=9 

ig=1 

ge), 1=1 500 

call atodtk(channel.ig.ival) 
110 zero=zero + 1val 

zero=zero/50().0/409.6 


write(* ,*) 
line=’Sampling finished. 
write(*,100) line 
memite( + .*) 
write(* ,*) 
line=’ Now please press the "-cal” button and hold it down’ 
write(*,100) line 
write(*,111) 
111 format(5x,’for 5 seconds. Hit return when ready’, 

A * to do this.’,\) 
read(*,104) ans 
write(*,*) 
write(* ,*) 
line=’Sampling simulated full scale speed and averaging’ 
write(*,100) line 


call setclock(0.01 ,0) 

full=0 

do 120 i=1,500 

call atodtk(channel,ig,ival) 
120 full=full + ival 

full=full/500.0/409.6 


write(* ,.*) 
line='Sampling finished.’ 


- 132- 





102 


104 


20 
30 


40 


> 


41 


write(*,100) line 
write(* ,*) 
write(* ,* ) 
write(*,102) zero, full 
FORMAT(5x,’Zero voltage: ’,F6.3,’ Full scale voltage: ’,F6.3) 

write(* ,* ) 

wnite(*, (A\)’)” Accept these? (y/n): ’ 
read(*,104) ans 
if(.not.((ans .eq. ’y’) .or. (ans .eq. ’Y’))) goto 99 
write(*,*) 
RETURN 

format(a) 
END 


Subroutine rpmcal 
common /rpm/ rvz, rvcal 
integer* 2 channel 
character ans* |] 


This subroutine is used to calibrate the Daytronics 
channel B, used to record propeller rpm data 


srite(* .*) 
svrite(* ,* ) 
write(*,10) 
format(10x,’ Welcome to the RPM (Channel B)’, 
’ calibration routine. ’) 
write(*,*) 
write(* ,*) 
write(*,20) 
format(15x,’A: Set Zero to read 000’) 
write(* 30) 
fonnat(15x,'B: Push cal button, adjust span to 1500’) 
write(*,*) 
write(*,*) 
pause ’ With propeller stopped, take time to do this now!’ 


write(* ,*) 
write(* ,*) 
write(* ,40) 
format(1x,'Hit any key while holding down channel B’, 
’ cal button for 5 seconds:’) 
read(*,41) ans 
format(a) 


channel=10 
rvcal=average(.05 channel) 


write(* ,*) 





50 


55 


> 


write(*,*) 
write(* ,50) 
format(1x,’ Please release cal button, and hit any’, 
’ key to take rpm zero: ’) 
read(*,41) ans 


rvz=average(0.05,channel) 


write(* ,*) 
write (*,55) rvz.rvcal 
format (2x,’Rvz = ’,f7.4,5x,’Rvcal = ’,f7.4) 
write (*,*) 
write(*,’ (a\)’) ’ Do you wish to recalibrate? (y/n): ’ 
read(* 41) ans 
if (ans .ne. ’n’) goto 5 
retum 
end 
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——— Version: 1.1 
ssoeee Date: 29 March 89 


------ Pupose: To take static balance readings at different points around 
the circle and store them for later use as a zero reference 
file. The data taken here represents the weight ( or W-B 
if done in water) effect AND the offsets of the balance. 


er Programmer: Dave Johnson 
------ Language: Microsoft Fortran 4.1 


------ Varibles: 
larray: row array in DMA filled by each MCBATOD 
call. 


zerot: array with averaged counts for all 
channels. This array will be read by 
the subroutine filter to convert to the 
preferred 6 channels and put into zeros(1). 


we @®@e eee 2 @® 2B ewe eee ew eB SSF 8 SS SFt es FSF SS SSS SE FSF FO FF 8 SH FF SSF SSS FF FS FFF FFF SF FFF FFF FFF eee 


integer* 2 1array(1000).bchan.echan,nsamp.nchan,n,nsect 
real zerot(32,9), angle, sum(9) 

character*11 fname, dcname 

character*3 tcase 

character*80 header 

character* 1 ans 

nsect = 32 

bchan = 0 

echan = 8 


rite (*,*) 
write (*,*) ’ ****** Welcome to DCZero! ********’ 
write (*,*) 


----- Input zero fuename 
write (*,*) 
write (*, ’ (A\) ’) ’ Type in name of zeroes file (11 char max): ’ 
read (*,’ (All), err= 1) fname 
open (unit = 1, fue = fname) 
write (*,*) 
write (*,*) ’ Type-in file header ...... (80 char max):’ 
read (*,’ (a80)’) header 
write (1,’ (a80)’) header 
write (1,2) 


—— full zerot(j.1) with zeros 


do 50 j=1,nsect 
do 50 1=1,nchan 
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zerot(j.i) = 0.0 


50 continue 
angle = 0.0 

C ----- Input Test case number 
write (*,*) 


write (*, ’ (A\) ’) ’ Type in test case (3 char. max): ’ 
read (*,’ (A3) ’) tcase 


C ------ Setclock set for 100Hz 


call mbopen 
call setclock(.01,1) 
write (*,*) 
write (*,’ (A\) ’) ’ Input Number of samples (e.g. 100): ’ 
read (*,6) nsamp 
write (*,*) 
20 pause ’ Press enter to start taking of zeroes ’ 


call mcbatod (bchan,echan,] ,nsamp.iarray ) 


c ----- calculate number of channels sampled 
nchan= (echan - bchan) + ] 
k=0 
n=n+] 

C ----- read in data to zerot(n.1) 


ae do 60 1=1, nchan 
zerot(n,i) = zerot(n,i) + 1array(i+(nchan*k)) 


60 continue 
k=k + ] 
if (k .EQ. nsamp) goto 71 
goto 55 
Wl continue 
C ----- Next, divide the values in zerot(n,i) to get averaged values 


do 80 1=1 ,nchan 
zerot(n,1) = (zerot(n,i))/nsamp 


80 continue 

C ----- write zerot(n,i) to file fname and to screen...... 
write (1,7) angle, (zerot(n,i) , i=1 ,nchan) 
write (*,*) 


write (*,2) 
write (*,7) angle, (zerot(n,i), i=1 ,nchan) 
write (*,*) 


angle = angle + 1)].25 





SA 


write (*,’ (a\)’) ’ Type dc counts output file name: 


-- Check to see if another mun 1s next .... 


if (n .eq. 32 ) then 
write (*,*) 
write (*,*) ’ N = 32 !, One full revolution has been done’ 
write (*,*) 
endif 
write(*,’ (a\)’) > Another run? (y/n):’ 
read (*,’ (al)’) ans 
if (ans .eq. ’y’) then 
goto 20 
endif 
close (unit=1) 


-- Next sum up all points/chan and put in dc ‘zero’ file .... 
do 90 1=1,9 


sum(i)=0.0 
continue 


do 96 i=1,9 


do 95 j=1,nsect 
Ssum(i)=sum(i) + zerot(,1) 
continue 
sum(i)=(Sum(1))/nsect 
continue 


-- Write to file ... 


write (*,*) 

read (*,° (al 1)’) dcname 

open (unit=2,fule=dcname,status= ‘new’ ) 
write (*,*) ° Input header for summed file ..... 
read (*,’ (a80)’) header 

write (2,’ (a80)’) header 

write (2,3) 

write (2, 4) tcase, (sum(1), 1=1,9) 

close (unit=2) 


-- FORMAT statements 


femmat ( Anple’ 3x, Ni’,5X,’N2’,5X, Y1’,5X, Y2’,5X,’RI’, 
Die (ou, OC) OX, AZ OM, XS ) 

format (‘tcase’,1x,’N1’,5X,’N2’,5X,’Y1’,5SX,’Y2’,5X,’RI’, 
ao eo, A OX, XZ OX, XS ) 

format (a3 ,9f7.1) 

format (14) 

format (f6.2,9F7.1) 


end 
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$DEBUG 


C ----- Version: 1.3 

C ----- Date: 27 April 89 

C ----- Programmer: Dave Johnson 

C ----- Language: Microsoft Fortran 4.1 


C ----- Purpose: To perform a fast fourier transform on coning 
C motion test data. 


C ----- Subroutines: 


C FAST: Performs the FFT and outputs the Fk’s 


complex x(320) 

integer npow,ndim,nreal nread 

real one, xreal,mag(320),freq(160).delta.rpm 
character*11] fname, outname,gname 


one = -1.0 
ndim = 320 


C ----- Input NPOW, the power of 2 that equals the # of data points .. 
write(*,*) 
write(*,’ (a\)’) ’ Input power, n, such that 24n = # data pts: ’ 
read (*, 5) npow 
nread = 2**npow 
write(* ,*) 
write(*,’ (a\)’) ’ Input the RPM data was taken at: ’ 
read(*, 6) rpm 
delta = 60/(32*rpm) 
write(*,*) delta 


¢ ---- Input Data file name ..... 
write (*,*) 
write (*,’ (a\)’) ’ Input name of data file: ’ 
read (*,’ (all)’) fname 


c ---- Read in data from file fname ..... 
open (unit=1 ,err=20,file=fname,status='old’) 
do 15 1=1,nread 
read (1,*,end=25) xreal 
x(i)=cmplx(xreal,0) 
15 continue 


goto 30 


20. ~write(*.*) ’ Error in opening file’ 
go to 100 
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25 write (*,*) ’ Error in reading file’ 
goto 100 
30 continue 


¢ ---- Call FFT subroutine.... 
call FAST(x,ndim,npow,one) 


c ---- Write FFT data to output file .... 
write(*,*) 
write(*,’ (a\)’) ’ Input FFT output fule name: ’ 
read (*,’ (all)’) outname 
open (unit=2,fule=outname,status=’new’) 
write(2,10) outname 
write(2,9) 
do 40 i=1 ,nread 
mag(i)=cabs(x(1)) 
wrnite(2,12) x(1),mag(i) 
40 continue 


¢ ---- Create data file for Grapher plot ...... 
write (*,*) 
write (*,’ (a\)’) ’ Input Grapher Data file name: ’ 
read (*,’ (al1)’) gname 
open (unit=3.flle=gname,status= ‘new’ ’) 
write (3,13) 0.,mag(1) 
do SO 1=2, (nread/2) 
freq(1)=(1-1 )/(nread* delta) 
write (3,13) freq(1),(2*(mag(1))) 
50 continue 
write (3,13) (freq(nread/2)+freq(2)), mag(nread/2 + 1) 


c ---- Fonnat statements .... 

5S format(13) 

6 format(f5.0) 

9 format(6x,’Real’,9x,’Imag.’,14x,’Magnitude’) 
10 format(2x,’ FFT Output File:’,A11) 

12 format(2x,f12.3,4x,f10.3, 9x, f12.3) 

13 format(2x ,f6.3,5x,f12.3) 


close(unit=1 ) 
close(unit=2) 
close(unit=3) 


100 continue 
end 
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c This is the fft routine. I’ll give you the meaning of ’ONE 


c later... 
SUBROUTINE FAST(X,NDIM,NPOW ,ONE) 


C 
C IMPLEMENT THE FAST FOURIER TRANSFORM 
C 


IMPLICIT REAL(A-H,O-Z) INTEGER(I-N) 
COMPLEX T,U,W,X(NDIM) 
DATA PI/3.1415926535897932384/ 
N1=2**NPOW 
N2=N1/2 
N3=N1-1 
J= | 
DO 300 I=1,N3 
IFU.GE.J) GO TO 100 
T= X(J) 
X(J)=X(1) 
X(1)=T 

100 K=N2 

200 IF(K.GE.J) GO TO 300 
J=J-K 
K=K/2 
GO TO 200 

300 J=J+K 
DO 500 L=1 .NPOW 
ei=2**L 
le 2=LE1/2 
ANG=PI/LE2 
U=CMPLX(1..0.) 
W=CMPLX(COS(ANG),ONE*SIN(ANG)) 
DO 500 J=1,LE2 
DO 400 I=J,N1,LE1 
IP=I+LE2 
a= X(IP)*U 
X(IP)=X(1)-T 

400 X(D= X(1I)+T 

500 U=U*W 
IF(ONE.EQ.1.)RETURN 
SCL= 1./REAL(N1) 
DO 600 I=1,N1 

600 X(I)=Xd)*SCL 
RETURN 
END 
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